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An  analytical  study  such  as  the  one  here  presented  is  of 
practical  value  only  if  it  can  be  readily  used  to  obtain 
uuaerical  results-  This  analysis  escapes  the  fate  of  being  a 
sere  nathesatical  exercise  because  Dr.  Joshua  B.  Greenspan,  of 
J  G  Engineering  Research  Associates,  Beltioore,  Maryland,  has 
applied  to  this  analysis  his  vast  experience  in  evaluating  the 
cocplicated  integrals  characteristic  of  radiation  problems  in 
cylindrical  coordinates.  These  quantitative  results  will  be 
presented  by  Dr.  Greenspon  in  a  companion  report:  "Axially 
Syiaetric  Green V  Functions  for  Cylinders."  Drs.  Alexander 
Silb’gnr  and  Swald  G.  Eichler  of  this  flrn  contributed  useful 
coc=»*ncs . 
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The  radiation  tape  dance  of  s  cylindrical  sound  source  of 
finite  length  can  be  expressed  as  the  sus  of  two  components: 

2  =  2  +  Z 
r  a 

where  Zp  ir  the  impedance  evaluated  bj  Keans  o '  the  generally 
used  technique  originated  by  Robey,  which  essvnes  that  tne  radi¬ 
ating  surface  is  bracketed  between  two  rigid  seed- infinite,  cy¬ 
lindrical  baffles.  2^.  is  associated  with  the  radisl  velocity 
distribution  a(t)  over  these  two  seed- infinite  cylindrical  sur¬ 
faces  and  is  therefore  in  the  n.ture  of  a  correction  factor  to 
Robey’s  inpedsnee.  a(z)  is  an  n known  function  which  satisfies 
a  non-homogeneous  Fredholm  integral  equation.  A  functional  J[a] 
is  contracted  which  is  stationary  and  proportional  to  for  the 
correct  solution  a(a): 

&J[a]/6a  =  0,  ZQ  *  J[a] 

Froo  this  vnriational  principle  a  value  of  is  calculated  by 
means  of  a  RayleJ gh-Ritz-type  procedure.  Finally,  the  far  field 
is  evaluated.  The  variational  principle  used  here  parallels 
the  Lerine-Scirvinger  principle  widely  used  to  obtain  scattering 
cross  sections.  Variational  solutions  are  presented  for  solid 
and  free- flooding  cylinders  for  axisynaetric  and  for  arbitrary 
velocity  distributions.  A  variational  solution  ic  given  for 
"squirters"  of  finite  wall  thickness,  but  It  is  ret  -ricted  to 
thin  walled  transducers.  In  an  Appendix,  non- variational  solu¬ 
tions  of  the  integral  equation  for  a(z)  are  presented  for 
"squirterr"  of  greater  vail  thickness. 

In  a  companion  report,  the  procedures  developed  here  ere  applied 
to  t’«  evaluation  of  the  self-  and  nutusl-radiaticn  impedances 
of  e.'-enects  in  an  array  of  coaxial,  free-flooding  axially  spaced 
ring  transducers .  Dr.  J.  S.  Greenspon,  cf  J  0  Engineering 
Ilssearch  Associates,  evaluated  the  inverse  Fourier  transforms 
required  for  these  solutions  and  obtained  quantitative  results 
which  he  will  present  in  s  separate  report. 


Systole'' 

(Alteras* ive  subscripts,  viz  u^,  are  used  to  condense  two  equations 

into  one,  the  upper  subscript  on  the  left  side  of  the  equation  being 
associated  with  the  upper  signs  and  subscripts  on  the  right  side  of 
the  equation,  and  vice  versa.) 

a  radius  of  cylindrical  source  (Fig.  2);  moan  radius  of  "squirter”  (Fig.  3) 
ai,aQ  inner  and  outer  radius  of  "squirter,”  respectively  (Fig.  3) 
c  sound  velocity  in  fluid  medium 

Gi,Go,gi,go,G^+,G1_,»'(z-z,)>  ana  H(s)  Green's  functions  anC  related  functions 
defined  in  table  1,  p.  1 

H  Hankel  function  of  the  first  kind,  of  order  a  (with  this  notation,  a 
n  massive  reactrace  is  negative) 

h  half  thickness  of  "squirter"  (Fig.  3) 

J  Seasel  function  of  order  m 

□ 

k  wave  number,  equal  to  <c/c 

kr  radial  wave  number,  equal  to  (k2-k^)^ 

axial  wave  number 

L  half  length  of  cylindrical  radiator  (Figs.  2  and  3) 

p  sound  pressure 

r,®,z  cylindrical  coordinates 

R,{?  spherical  coordinates 

U  radial  velocity  amplitude  of  cylindrical  source  (Fig.  3) 

u(z)  radial  velocity  distribution  on  cylindrical  surface  (r=a),  positive 
outward  (Fig.  3) 

Z  radiation  dope dance  of  cylindrical  source  in  units  of  force/velocity 

equal  to  (Zj+Z^)  (^iS-  2) 

7  ,  radiation  impedance  obtained  from  Robey's  model  (Pig.  lb)  for  which 
a(z)  =  0 

Z  correction  factor  associate  :  with  vsjocity  distribution  o(z)  and  to  be 
-dded  to  Z  (Fig.  2) 

a( z)  radial  velocity  in  the  regions  ,z|  >  L  (Fig.  2)  normalized  to  velocity 
amplitude  U  of  cylindrical  surface 

»  Poisson's  ratio  of  transducer  material 

0  density  of  fluid  medium 

V  velocity  potential  (outward  velocity  18  -5$/dr);  subscripts  "i"  and  "o" 
refer  respectively  to  regions  r  <  a  and  r  >  a 

u>  circular  frequency  (harmonic  time  dependence  factor  exp(-icot)  which  multi¬ 
plies  the  velocities  and  the  potentials,  has  been  suppressed  throughout 
this  report] 

*OtUor  symbols  are  defined  in  the  text  . 
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IV. 1 

m 

All  configurations  with 
axlsynasetric  velocity 
distribution 

General  form  of 
axisynzaetrlc 

Green's  function 

6^,8  ,z-z*) 

■ 

■ 

Free- flooding  cylinder  with 
axJsyatLStrlc  velocity 
distribution 

8o0*,6,z-z') 

IV.  6 

i  r  >  a 

l 

I 

j 

All  configurations  with 
axisysaaetric  velocity  dis¬ 
tribution  symefitrical 
about  z  =  0 

Even  component  of 

Go 

VH2.2 

r  <  a 

Free- flooding  cylinder  with 
axisysmftric  velocity  dis- 
tributi'«i  sysssotrical  about 
z  =  C 

r-~ — 

i  Even  component 
of  Gi 
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2-Z  } 

IV. 3 
and  b 

r  <  a 
z  *>  L 

Solid  cylinder  with 
axisysssetric  velocity 
distribution 

Green's  function 
whose  normal  de¬ 
rivative  d/dz 
vanishes  cu  end 
cap  z  =  L 

°'-(Ul’> 

IV.  3 
and  b 

_ _ 

■ 

Green's  function 
whose  normal  de¬ 
rivative  d/dz  van¬ 
ishes  on  end  cap 
z  =  -L 

°oa(r‘'j 

z-z  ) 

■ 

Cylinder  with  arbitrary 
velocity  distribution  ex¬ 
pressed  as  a  Fourier  aeries 
in  ©  (Bq.  X.1) 

General  fora  of 
non-axisymaetric 
Green's  function 

r!'5 

r  <  a  | 

J 

r(z-z') 

r  *  a 

Solid  cylinder 

Linear  combination 
of  Green's 
functions 

VUI.ba 

and 

van. 5 

r  *  a 

Free-flooding  cylinder  of 
vanishing  wall  thickness 

IX.lba 

r  c  a 

"Squirter" 

n(z) 

IV. 8 

Solid  cylinder 

Integral  of  Green's 
function  war 
radiating  surface 

VIII. bb 

Free- flooding  cylinder  of 
vanishing  wall  thickness 
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I. 


In  this  report  a  variational  technique  is  used  to  derive  expressions  for  the 
radiation  loading  of  radially  pulsating  cylinders  of  finite  length.  Bid  effects 
are  accounted  for,  no  restrictions  being  placed  on  the  circulation  of  the  acous¬ 
tic  fluid  around  the  edges  of  the  cylinder.  The  far  field  potentials,  on  the  axis 
of  the  cylindrical  radiator  and  for  other  bearings  are  also  given.  The  analysis 
is  performed  for  (l)  a  solid  cylinder,  (2)  an  open-ended,  free-flooding  cylinder 
of  vanishing  vail  thickness,  and  (3)  a  "squlrter"  of  ora 11,  but  non-vanishing 
thickness-to- radius  ratio.  The  final  section  presents  an  extension  of  the  t  al- 
ysis  to  cylindrical  radiators  embodying  an  arbitrary,  non-axisymaetric  velocity 
distribution.  A  non-variatianal  technique  presented  in  Appendix  B  extends  this 
study  to  "squirteio"  01  1'rger  vail  thickness. 

With  the  present  approach,  the  need  for  machine  calculations  has  been  con¬ 
fined  to  the  evaluation  of  Green's  functions  in  the  fora  of  Robey's  integrals. 

As  mentioned  in  the  acknowledgment ,  these  integrals  are  being  evaluated  by  Dr. 
Greenspon,  J  G  Engineering  Research  Associates.  The  technique  developed  in  this 
study  has  been  extended  to  the  evaluation  of  the  mutual  radiation  impedances  be¬ 
tween  elements  of  an  array  of  free- flooding  ring  transducers.0,  Humeri cal  results 
for  this  configuration  are  also  being  obtained  by  Dr.  Greenspon  and  will  be  in¬ 
cluded  in  his  report:  "Axially  Symmetric  Green's  Functions  for  Cylinders.” 


II.  A  Review  of  Published  Analytical  Approaches  to  the  Finite  Cylindrical 
Radiator 


The  fluid  potential  $  generated  by  a  sound  radiator  is  given  by  the  familiar 


Helmholtz  integral  equation 


«(5)  -  j  [G(R,R')  *(g')]  ds' 

t*  t 


(rt.i) 


*73xis  analysis  will  be  presented  in  a  report  to  be  published  in  March  196k, 
"Mutual  Radiation  Impedance  fo-  Spaced,  Coaxial,  Free-Flooding  Ring  Transducers," 
CAA  Report  U-178-^9,  Contract  Konr-2739- 
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where  8'  is  the  radiating  surface,  and  n'  is  the  outward  noraal  to  the  surface  of 
integration.  The  first  terra  in  the  integrand  can  be  readily  evaluated,  if  the 
noraal  velocity  u(5'),  of  the  radiating  surface,  which  equals  -dtySa',  is  known. 

The  second  tens  in  the  integrand  involves  the  unknown  potential  on  the  radiating 
surface,  4(5').  We  Bust  therefore,  in  general,  solve  an  integral  equation  to  ob¬ 
tain  this  potential.  In  the  last  two  years,  the  general  availability  of  large, 
digital  computers  has  made  it  practical  to  use  a  finite-difference  method  to  ob¬ 
tain  numerical  solutions  of  the  Helmholtz  integral  equation  for  the  finite  cylin- 
drical  radiator  (Fig.  la  and  the  Table  on  p.3).  The  drawback  of  this  approach 
is  that  the  large  computational  effort  involved  mist  be  repeated  for  every  combi¬ 
nation  of  length-to-radius  ratio,  of  ka,  tad  of  surface  velocity  distribution. 

Another  successful  approach  to  the  finite  cylinder  problem,  which  circumvents 
the  Helmholtz  integral  equation,  uses  an  expansion  of  the  potential  in  spherical 
wave  harggri.es.  For  this  approach  the  volume  of  calculations  is  less  than  for 
the  finite-difference  approach  described  above.  Thus,  approximate  results  were 
obtained  in  ref.  4  without  the  help  of  electronic  computers,  by  confining  the 
series  expansion  to  only  a  few  terms.  For  practical  applications,  this  method 
also  requires  computer  facilities. 

An  approximate  method  which,  historically,  precedes  the  approaches  described 
above,  consists  in  constructing  a  Green's  function  whose  derivative  dO/Sa'  van¬ 
ishes  mi  the  infinite  cylindrical  surface  r=a.  If  we  now  prolong  the  cylindri¬ 
cal  raaiator  by  two  seal-infinite  rigid  cylindrical  baffles  of  the  same  diameter, 
the  surface  integral  in  Eq.  XX. 1  is  confined  to  the  cylindrical  surface  (Fig.  lb) . 

Over  this  surface,  the  second  term  in  the  integrand,  which  involves  the  unknown 
potent 'al,  has  been  eliminated  by  our  choice  of  the  Green's  function.  We  can 
therefore  obtain  an  approximate  expression  for  the  potential  without  having  to 
solve  an  integral  equztiza: 

4(5)  »  -  J u(r')g(r,r')  as'  (ix..O 

With  this  approach  laird  and  Cohen"’  derived  an  analytical  expression  for  the  far 
field  potential,  the  integrals  being  evaluated  by  the  method  of  stationary  phase. 
These  integrals,  which  for  the  axisyametric  velocity  distribution  are  known  as 
Robey's  integrals,  must  unfortunately  be  evaluated  mnericolly  if  the  potential 
on  or  near  the  radiating  surface  is  required.^  Greenspan  has  simplified  the 
technique  ror  performing  this  integration. 1  He  and  Shorten  also  evaluated  these 
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integrals  for  noa-exisynaetrie  velocity  distributions .  The  drawback  of  Robey's 
mathematical  aodel  is  that  it  does  not  permit  circulation  of  the  fluid  around  the 
edges  of  the  transducer,  because  of  the  assumption  of  two  semi-infinite  cylindri¬ 
cal  baffles.  Neither  does  this  aodel  lend  itself  to  the  evaluation  of  axially 
vibrating  solid  cylinders,  or  of  the  free-flooding  open-ended  cylinders  known  in 
acoustical  vernacular  as  "aquirters.®  Robey  originally  approximated  the  radiation 
loading  or  such  free-flooding  transducers  by  assuming  that  the  fluid  colunn  inside 
the  "squirter"  is  terminated  by  pressure-release  piston a.^  He  then  refined  his 
analysis  by  assuming  the  terminal  impedance  to  be  that  of  £  piston  in  an  infinite 
plane  baffle10  (Pig.  lc). 

In  summary,  existing  analyses  use  either  a  large  computational  effort  which 
must  bo  repeated  for  every  particular  combination  of  sound  source  parameters,  or 
sa  elegant  approximate  technique  which,  however,  does  not  account  for  the  circu¬ 
lation  of  fluid  around  the  two  extremiti-..  of  the  cylinder. 


m.  Description  of  the  Present  Approach 

She  present  approach  makes  use  of  a  Green's  function  similar  to  Robey* 3, 
tijs  eliminating  from  tbs  deluholtz  equation  the  tern  containing  the  potential  on 
the  cylindrical  surface  r=a.  However,  instead  of  assuming  the  source  to  be 
bracketed  by  rigid  beffles,  the  potential  is  expressed  in  tenas  of  an  unknown 
radial  velocity  distribution,  a(z),  over  the  two  seal'- infinite  cylindrical  bound¬ 
aries  prolonging  the  sound  source  (Fig.  8).  The  potential  in  the  cylindrical 
column  in  the  rsgion  r  <  a  is  then  formulated  with  the  help  of  a  suitable  Green's 
function  whose  normal  derivative  vanishes  on  the  boundary  r*a.  The  potential  in 
this  cylindrical  regies  in  also  expressed  in  terms  or  the  unknown  velocity  distri¬ 
bution  a(z).  By  requiring  continuity  of  these  two  potentials  across  the  cylindri¬ 
cal  boundary  r= a,  |z ]  >  L,  a»:  integral  equation  for  a(z)  is  obtained.  If  we  coo¬ 
ps ’•e  this  formulation  to  the  free-space  Green's  function  formulation  in  ref.  2  and 
3,  we  see  that  t»  unknown  function  a(s)  and  the  surface  r*a,  Jzj  >  L  take, 
respectively,  t is  place  of  $(H/)  and  of  the  radiator  surface.  One  advantage  of 
the  present  approach  is  that  a  gives  error  in  the  expression  ter  a(s)  esn  be  ex¬ 
pected  to  result  in  a  smaller  error  in  the  n.  Uation  tspedanee  than  would  result 
from  a  siailer  eiror  in  in  ref.  2  and  3.  The  principal  advantage,  however, 

is  that  this  approach  lends  itself  to  an  aoprcecisate  varieties!  1  solution  both  of 


the  solid  and  open-ended  finite  cylinder.* 

She  radiation  impedance  Z  can  be  written  as  the  sum  of  the  impedance  ob¬ 
tained  by  setting  a(z)  =  0,  and  of  an  impedance  ZQ  associated  with  a(z),  the 
unknown  velocity  distribution  in  the  region  (2  |  >  L: 


Zr+Sa 


(2n.i) 


Z  is  thus  in  the  nature  of  a  correct ioa  factor  to  ispedances  computed  by  Robey, 

a  7  8  9 

dreenspcn,  and  She roan  fron  Robey's  oathesstical  model.  By  virtue  of  the 
variational  principle  to  be  derived  in  Section  V,  for  the  correct  solution  of  the 
integral  equation  a(z),  is  proportional  to  a  functional  J|ct*  which  is  station¬ 
ary  with  respect  to  first  order  variations  of  a(s) : 


-la] 


isisi  =  0 

60e 


(Xtt-2) 


Purthernore ,  J[a]  depends  on  the  functional  forn  of  a(z)  but  not  on  its  anplitude. 
The  technique  for  confuting  is  similar  to  the  Rayleigh-Rltz  technique  for 
evaluating  the  natural  frequency. 

This  approach  jerellela  the  use  of  the  Levine -S cbuinge r  variational  principle 

fo?  bcoitering  cross  sections,  which  has  beer  applied  to  a  large  number  of  dif- 
12 

fraction  prowlers.  The  equivalent  principle  for  radiation  Impedances  is  proved 

11 

in  its  general  fora,  U3irg  free-space  Greek's  functions,  by  Morse  and  Feshbacfc. 

These  authors  do  not,  however,  use  it  to  solve  any  particular  problem.  Apparent¬ 
ly 

ly,  only  Store r  applied  this  principle  to  a  specific  problem,  viz.  tn»  effect 

of  a  finite  circular  baffle  on  the  radiation  loading  of  a  coaxial  antenna.  In 
195 y,  Professor  Storer,  of  Harvard  University,  suggested  to  the  author  of  this 
report  t!,at  the  axlsyBnetrically  vibrating  cylinder  of  finite  length  could  also 
be  analyzed  in  this  fashion.  Consequently  a  rather  sketchy  variational  solution 


*A  rigorous  Wiener-Hop?  type  solution  of  the  Integral  equation  is  oossible  for 
seal-infinite  cylindrical  radiator  problems  formulated  in  this  fashion,  accord¬ 
ing  to  Levine  and  Schwinger.11  Those  authors  mention  this  formulation  as  an  al¬ 
ternative  to  the  one  they  actually  used  in  their  analysis  of  sound  radiation  from 
a  semi- infinite  pipe.  Levine  ,l26  extended  this  study  to  pipes  of  arbitrary  cross 
section.  One  of  his  approximations  for  the  reflection  coefficient  is  obtained 
from  a  variational  solution  of  an  integral  equation  (his  "variational  principle  A") 
a pp liable  over  the  seal-infinite  cylindrical  surface  extending  the  pipe,  and  is 
therefore  of  the  fora  of  the  Integral  equation  used  in  this  report. 
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of  the  solid  cylindrical  source  vas  presented  In  an  Internal  aesorandum  of  the 

15 

Harvard  Acoustics  Research  laboratory.  A  detailed  analysis  vas  not  carried  out 
because  numerical  results  depended  on  the  evaluation  of  Robey's  integral,  which  at. 
that  tine  was  not  available.  Since,  as  mentioned  earlier,  such  integrals  can  now 
be  readily  evaluated,  and  since  Dr.  Graeaspos  hindly  egr*ed  to  apply  his  experi- 
ience  in  this  type  of  calculation  to  tha  problem  at  hand,  it  is  now  worthwhile  to 
use  the  variational  formulation  to  obtain  a  solution  to  the  finite  cylinder 
problem. 


IV.  In teg 


auatlon  Formulation  of  the  Solid  Cylinder  Problem 


The  infinite  region  surrounding  the  cylindrical  radiator  is  subdivided  into 
three  regions  (Fig.  2):  an  outer  region,  r  >  a,  identified  by  subscript  o;  and  two 
semi-infinite  inner  regions,  r  <  a,  one  corresponding  z  >  L  and  identified  by  the 
subscript  i+,  and  a  second  inner  region  corresponding  to  values  of  z  <  L,  identi¬ 
fied  bv  the  subscript  i-.  The  Green's  function  for  the  outer  region  satisfying 
the  condition  30/dr#*O  for  r*=c,  was  constructed  by  Robey 


.  f  H  (k  r) 

Go(r,a,z-z')  .  -jr  j  expli^fz-z')]  a**  (IV 

•39 

2  2  4 

She  evaluation  of  G  is  the  subject  of  references  6  and  7-  Since  k*(k -*.)*,  the 
o  r  z 

functions  or  x  in  tbs  integrand  ore  even  in  Hence,  only  the  real  component 
of  the  exponential  function,  costh£(z-z*)J  contributes  to  the  integral.  The  seme 
comment  applies  to  the  Green's  function  for  the  infinite  cylindrical  region  r  <-  a: 


(iV.l) 


(r,o,z-z*)  - - §-  f  —  exptik  (?-z*)3 

1  kao  £  k^Oy*)  s 


(IV.2) 


This  function  is  derived  and  evaluated  in  Appendix  A. 

In  the  csss  of  tha  solid  cylindrical  radiator  it  is  convenient  to  coatoins  two 
Green's  functions  of  the  fora  of  Eq.  IV.2  so  that  the  normal  derivative  of  the  re¬ 
sultant  Green's  function  vanishes  over  the  end  capo  of  the  cylinder,  i.e.,  in  the 
two  circular  regions  r  <  a,  i  «  il.  This  will  eHjdrate  the  potential  tern  fro® 
the  Bslsholtz  integral  owr  the  end  caps,  as  well  as  aver  the  cylindrical  surface. 


V 


Such  Green's  functions  are  readily  constructed  by  introducing  isage  sources: 


Gi+{r,8,S~Z*^  =  +  Gi(r»a»z42'-2L) 

G^tr.SjZ-z')  =  G^rj&jZ-z')  +  Gj(r,a,z+z'4  2L) 

These  Green's  functions  can  be  written  core  concisely  as* 


Git<r,a ,*-*') 


(IV.3) 


(IV-M 


We  can  now  write  the  potentials  in  these  three  regions  by  caking  use  of  the  codi¬ 
fied  Helnholtz  integral  in  Eq.  H.2.  If  ve  assuas  that  the  end  caps  ere  rigid, 
the  surface  integral  reduces  to  the  cylindrical  surface:** 


$o(r,z)  =  2xa  J  u(z,)Go(r,a,z-z')  dz' 


(IV.' 


$i+(r,z)  =  +5sa^  u(z,)tG^(r,&,z--£")  +  Gj(r,a,z4-z‘  +  2L)]  dz' 


(IV. 5b) 


The  tine  dependence  of  u(z*)  and  of  the  potentials  is  hamonic-  The  tine-dependent 
function  exp(-ia)t)  has  beer,  emitted,  for  the  soke  of  brevity  throughout  this 
report.  These  integrals  are  of  opposite  sign,  because  dfydn'  in  the  Helnholtz 
equation  equals  -di^/dr^u  in  the  outer  region,  and  +£a./dr=  -u  in  the  inner  region. 
With  this  sign  convention,  the  sound  pressure  equals  p|.  If  the  velocity  distri¬ 
bution  of  the  radiator  is  oytEetriual  about  the  plane  2=0,  the  two  inner  regions 
will  have  identical  potentials  end  we  need  concern  ourselves  with  only  one  inner 
region  which  we  will  designate  by  the  subs'—ipt  i.  Pot  the  case  of  a  syazaetrlcal 
velocity  distribution,  only  the  part  of  the  Green’s  function  which  is  symetrical 


♦Here ,  end  elsewhere  in  this  report,  alternative  subscripts  and  signs  have  been 
used ,  for  the  sake  of  brevity,  to  condense  two  equations  into  one,  the  upper  sub¬ 
script  on  the  left  side  of  the  equation  being  associated  with  the  upper  sign  on 
the  right  side  of  the  equation,  and  vice  versa. 

**12  section  X  expressions  are  given  for  an  arbitrary  velocity  distribution  over 
the  radiating  surface.  The  potential  contributed  by  vibrating  end  caps  is  given 
in  Bq.  X. 1  and  5. 
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about  z'«0  cariribites  to  tbs  potential.  We  sill  denote  this  even  coaponent  or 


ao  V 


,  „  if  Ho<V> 

Bo(w-z)«;7rJ 


cos  k_z  ccs  k  a'  dk 
2  2  2 


(IV.6) 


We  further  specialise  the  problem  by  assusing  that  the  radial  velocity  of 
the  sound  source  Is  constant  and  equal  to  0.  0(2')  is  therefore  a  known  function 

for  |z {  *  L.  The  velocity  distributico  along  the  cylindrical  boundaries 
2  >  L,  r*a.  Is  an  unknown  function,  say  tfcr(i).  She  Integrals  for  the  potentials 


now  beouee 


^(r.z)  =  -2na  Tlj  afz'JtO^r.ajZ-z')  +  Gi(r,a ,a*s '-2L))  dz' 
L 


(I?.7a) 


*o(r,2)  =  fcsa  0{  j’  gQ(r,a,2-2')  dz'  ♦  Ja(s')go(r  ,8,2-2')  dz']  (lV.7b) 

0  L 

For  the  sake  of  brevity  we  will  fro*  new  on  express  the  known  coapoaent  of  the 
potential  5q,  evaluated  on  the  cylindrical  surface  r=a,  as  a  f  nction,  say  H{z) , 


rather  than  as  an  integral 


r 

H(  i)  «  J  g0(a,a. 


z-z')  d z‘ 


(IV. 8) 


We  can  now  construct  the  integral  equation  which  the  unknown  function  ct(z) 
cost  satisfy.  This  integral  equation  is  derived  froo  the  require*  nt  that  the  po¬ 
tentials  be  continuous  across  the  cylinder  boundary  r=a,  t  >  L: 


»0(a,z)  -  it(a.*)  0  0,  for  z  >  L 


(IV.9) 


Wnen  we  substitute  Eqs.  17.5,  this  continuity  condition  takes  the  fora  of  s  non- 

u _ - a» —  — «  -  4 ^ a  4—-  fSa  f{ra4> 

IMBUQCUCVUO  riqw.-,rin  7JMV4%4«  </a  -- —  iii3w  ftlUv* 


r(2-2')  a(z')  ds'  •>  -3(2),  for  2  >  L 


(XV. 10) 


where  the  kernel  01'  this  equation  is  given  by 

r{z-s')  =  go(o,aj2-a/)  +  |  G^a.a.z-z')  +  |  Gi(a,a,z+z/-2L)  (lV.il) 

Wa  will  now  show  that  the  function  a( z*)  which  satisfies  this  integral  equation 
gives  a  stationary  value  for  the  radiation  impedance,  with  respect  to 
variations  6a. 


V.  Derivation  of  the  Variational  Principle  for  the  Radiation  Impedance 

We  note  for  future  use  that  the  reailtant  radiation  impedance  on  the  cylinder 
is  obtained  by  integrating  the  pressure 

p(n,z)  =  pfc0(a,z) 


=  -hop$o{a,z) 

over  the  surface  of  the  cylinder: 


(V.l) 


z . .  asp  |  a* 

0 

J*  * 

=  -  i(«sa)2®p  J  [H(z)  +  J  o(z')eo(^,&,z-z')  dz *]  dz  (V.2) 

0  L 

The  radiation  inpedance  is  thus  clearly  the  sum  of  two  component  impedances,  as 
indicated  in  Eq.  III.l.  The  impedance  computed  from  Robey 1 3  model  is 


Zf  =  -(baa)  hup  j  H(z)  dz 


t 


1  L 

r  r 


=  -(baa)*"  imp  !  g  (a,a,z-z')  at'  dz 


(V-3) 


Like  H(z),  Is  a  known  quantity  since  it  does  not  involve  the  unknown  function 
a(z).  The  correction  tero  in  Eq.  ZZI.l,  which  esbodiesi  the  contribution  of  the 


9 


flow  across  the  cylindrical  boundaries  prolonging  tne  source  is 


/«/« 

0  L 


zn  =  -C1**®)  imp  [  0(2')  g  (a,a,j»-2/)  dz']  dz 


(V.4) 


Since  tbe  Green’s  function  is  syaectrical  in  2  and  z‘  the  order  of  integration  in 
Eq.  Y.4  can  be  inverted: 


Za  =  -(Ms*)  imp 


^  t/*e0(a»a>?-2/) 


]  a{z)  dz 


H(z)  a(z)  dz  (VO) 

A  functional  J(a]  will  be  defined  below,  Eos.  V.9  and  ID.  Ibr  future  reference 
we  note  that  for  the  correct  function  a(z),  i.c. ,  for  the  function  which  satis¬ 
fies  the  integral  equation,  Eq.  IV.  10,  this  functional  con  also  be  written  as 

CD 

j[a]  ■  -f H(z)  a(2)  dz,  for  a(z)  solution  of  Bq.  IV.10.  (V.6) 

L 

Comparing  this  with  Eq.  V.5,  we  can  relate  this  functional  and  the  iapedance  ZQ: 

p 

Za  =  imp(!*«a)  J(a],  for  a(z)  solution  of  Eq.  IV. 10.  (V.7) 

We  will  new  prove  that  J[aJ  is  stationary  with  respect  to  first  order  varia¬ 
tions  6a  about  the  correct  function  cr(z).  For  this  purpose,  we  relate  J[«]  to 
the  integral  equation,  Eq.  IV. 10.  Ve  nultiply  both  sides  of  this  equation  by 
a(z)  end  integrate  with  respect  to  z  over  tbe  region  z  >  I>.  We  then  divide  both 
sides  of  the  equation  thus  obtained  by 

if  H(z)  <s(z)  dz]2 
L 


=  -(4aa)  imp 


I 


Our  original  integral  equation  now  takes  Ine  forn 


-1 


a(z)r(z-z')  a(z')  dz*  dz 


[  ■  H(z)a(z)  dzf 


CD 

J  H(z)a(z 


Hie  functional  J[a]  is  defined  as  the  reciprocal  of  the  left  side  of  this  equation: 


Jia]  = 


O 

l[a]  =  J  H(z)  a(z 


(V.lOa) 


®  CO 

=  j  Ja(z)  r(z-z')  a( z‘)  dz'  dz 


(V.lOb) 


For  a  function  a(z)  vhich  satifi'ies  the  integral  equation,  and  therefore  the 
equality  in  Eq.  V.8,  the  reciprocal  of  the  right  side  of  Eq.  V.8  is  also  equal  to 
the  J[a],  as  already  indicated  in  Eq.  V.6.  If  the  functional  derined  in  Eq.  V.9  is 
indeed  stationary  vith  respect  to  snail  variations  of  the  function  0(7),  then,  by 
definition, the  increnent  6J[a]  associated  with  an  in ere sent  6a  is  zero: 


■»{ 


H(z)  6a(z)  dz 


_  AjISi  .  r(z-z')[a(z)6a(z'Mz')6a(z)l  dz'  dz  =  0  (V.ll) 

r.'ivo  the  Green's  functions  in  Eq.  IV. 11,  r(z-z')  is  synsctrical  with  respect  to 
z  and  z'.  We  can  therefore  i’-vert  the  order  of  intention  in  the  forcer  of  the 
two  terns  of  the  integrand  of  the  double  integral  in  Eq.  V.ll.  The  double  inte¬ 
gral  can  thus  be  condensed  to 


'  r  F  ? 

I  =  2  i  r(z-z')  ct(z)ba(z )  dz' 
'  J  *  J 
L  l 


(V.12) 
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I 

I 

I 

From  the  integral  equation,  Eq.  IV.  10,  ve  see  that  the  Integral  over  a'  equals 
-H(z)  for  the  correct  function  '3(2).  The  double  integral  can  thus  finally  be 
written  as 


jj“  -2J  H(s)  6a(z)  da  (V.13) 

L 

When  ve  substitute  Eq.  V.13  in  place  of  the  double  integral  in  Eq.  V.ll  and  mul¬ 
tiply  the  terms  of  this  equation  by  the  ratio  B^aJ/SAfci],  ve  obtain 


*Bla]  +  Ata]} 


a)  6a(z)  dz  «  0 


(V.lA) 


Since  the  integral  is  not  identically  zero  the  sum  in  brackets,  vhich  multiplies 
this  integral,  must  vanish,  i.e., 

B[a]  =»  -A[a]  (V.15) 

When  de  substitute  the  definitions  of  these  tvo  functionals,  Eqs.  V.10,  this 
bectses: 


CD  06  0» 

fU)  r(z-z')  a(z')  dz'  d z  *  -  f H(z)  a(z)  dz  (V.l£) 

j  J  J 

11  1 

This  equation  is  obviously  satisfied  if  cx(z')  satisfies  our  original  integral 
equation,  Eq.  IY.10.  We  have  thus  sbovn  that  the  functional  J[ot].  as  defined  in 
Eq.  V.y,  does  indeed  take  on  a  stationary  value  for  the  correct  value  of  the 
function  a(z) .  Since  the  functional  J[ct]  is  stationary  with  respect  to  the  cor¬ 
rect  function,  the  error  in  J[ctl  is  of  a  hi^er  order  than  the  error  in  a{s). 

We  shall  now  illustrate  the  evaluation  of  the  radiation  impedance  by  means 
of  the  variational  principle  Just  derived. 
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VI.  Evaluation  .f  the  Radiation  Impedance  from  the  Variational  Pilnciple 

We  first  proceed  to  select  the  simplest  trial  function  a(z)  which  yields  e 
far  field  potential  4^  in  the  desired  fora  of  a  spherically  spreading  wave, 

4j(a,s)  =  ~j-  exp(ikz)  ,  for  large  !s|  (VI.  1) 

If  the  cylindrical  boundary  r=a,  jz  |  >  L  were  in  the  fora  of  a  rigid  pipe,  the 
far  field  potential  would  only  decay  as  a  result  of  viscous  losses,  as  embodied  in 
an  ice  ginary  component  of  the  wave  number.  The  far  field  potential  in  a  rigid 
pipe  can  therefore  only  decay  exponentially.  The  desired  spherical  spreading 
loss,  Eq.  VI. 1,  must  therefore  be  the  result  of  energy  flow  across  the  cylindrical 
boundary  associated  with  the  velocity  a(z).  The  rate  of  energy  outflow,  per  unit 
axial  distance,  along  the  "pipe"  is 

=  akapcU  ^(ajzjjjajz)  j,  (VI.?) 

This  isust  balance  the  decrease,  per  unit  axial  distance,  of  *-hs  acoustic  energy- 
propagating  down  the  pipe: 


In  the  far  field,  and  for  the  values  of  kt  characteristic  of  transducers , 
can  be  set  equal  to  ^(a^).*  We  can  thus  replace  the  integral  in  Eq.  VI. 3  by 
a^ii(a,s)p/2.  When  we  substitute  Eq.  VI.  1  for  9^  in  Rqs.  VI.2  ncd  3,  ve  ssn  sc-lvv 
for  ja(z)  ,  : 

■a(z)|  =  »  for  large  jz!  (VIA) 

Is! 

Wc  thus  conclude,  that  in  the  far  field,  a(z)  must  be  of  the  fora 

a(z)  «  f  f0».  iargC  jZ  (VI.5) 

'  z 

We  will  now  illustrate  the  use  of  the  variational  principle  by  eelacting  t.he 

•Even  for  large  ka ,  the  value  of  4^  averaged  over  the  cylindrical  cross  se'tion 
is  equal  to  a  constant  tines  4^(o,z).  The  functional  relation  derived  in  Eq.  VI. b 
therefore  still  holds,  but  the  constant  in  this  equation  will  not  equal  A. 
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siaplest  trial  function  vnich  satisfies  Eq.  VI. 5  and  which  can  also  account  for  a 
no  re  rapidly  decaying  rear  field.  Such  a  fimction  requires  the  use  of  at  least 
two  unknown  coefficients,  and  Xgt 


X  x 

<*(z)  =  (~§  +  -f )  exp(lkz) 


(VI.6) 


Since  x1  and  Xg  are  generally  coeqilex  quantities,  this  expression  allows  for  a 
phase  shift  between  the  two  ccaponents  of  the  potentials.  By  virtue  of  the 
variational  principle,  the  best  values  of  the  unknown  coefficients  are  those  which 
give  a  stationary  value  to  the  functional  Jte]: 


=  0 

ajtaj/axg  =  0 


(VI.?) 


When  we  substitute  Eq.  V.9  for  the  functional  J[a]  we  can  write  these  equations, 
after  some  manipulation,  as 


a*felA[01_2§feljIal,o 


(VI. 8) 


etc. 

When  we  ccebine  the  assumed  trial  function,  Eq.  VI.6,  with  the  definitions  of  the 
functionals  A[a)  and  E{a],  Bqs.  V.10,  these  two  functionals  are  found  to  be, 
respectively,  linear  and  quadratic  in  and  Xg, 

A(a]  =  +  agXg,  dA/dXg  =  ag 

B[a)  =  b^  +  bgXg  +  2*^  x1xg 

dB/c =  abjX,*  oB/^Xg  =  SbgXg+  36^  x%  (VI. 9) 

where  coefficients  a^  ag,  b^  bg  and  b^  are  known,  complex  quantities.  When  ex¬ 
pressions  VI. 9  are  substituted  in  Eqs.  VI. 8  a  set  of  two  simultaneous  linear 
equations  is  obtained: 

2(a^  -  bx  JEaDx.  +  2(a ^g-bjg  Jfr])xg  =  0 

2^aia2  *  b12  J(alK+  2^82'b2  ”  0  (VI.  10) 


lb 


8 luce  this  set  of  equations  is  homogeneous ,  its  coefficient  matrix  must  vanish: 


1 

I 

7 

i 


] 

I 

i 


1 

I 

! 

i 


2{^-Vta])  2(a1VVrIaJ) 

^Vsfb^Jla] )  2(82-b2J^^ 


(VI.  U) 


Expanding  this  matrix  we  only  obtain  terms  proportional  to  J  [a]  and  Jfa].  The 
equations  can  therefore  be  solved  for  J[a]: 


J[a] 


•h  *  4i  -  fVjS 

v?  -  4 


(VI. 12) 


The  remarkable  feature  of  this  result  is  that  the  functional  ia  independent  of 
ana  x^.  Like  the  coefficients  a^,  ag,  b^,  bg  and  b^  the  functional  in  Eq.  VI. i? 
is  coaplex.  We  see,  by  referring  to  Bq.  V.7,  that  the  imaginary  component  of 
J[a]  embodies  the  radiation  resistance,  and  its  real  component,  the  reactance. 

If  we  are  merely  interested  in  the  radiation  loading  we  need  not  evaluate  the 
unknown  coefficients  and  .  If,  however,  we  wish  to  compute  the  far 

field  potentials,  we  must  substitute  the  expression  for  a(s)  in  Eqs.  IV. 8,  Bnd 
therefore  require  the  values  of  x^  and  Xg.  The  ratio  of  these  two  coefficients  is 
obtained  from  either  of  the  two  homogeneous  equations,  Eq.  VI. 10 


*2  (-a?  ♦  bjJlal) 

X1  3  Vs  •  bi7^ 


(VI. 13) 


where  the  value  of  J(a]  is  known  from  Sq.  VI.  12.  The  coefficient  x^  is  obtained 
by  substituting  the  ratios  Xg/x^,  uq.  VI. 13  in  Sq.  VI.6,  which  is  then  substituted 
for  ct(z')  in  the  integral  equation,  Eq.  IV.10.  Unless  the  functional  dependence 
or  the  trial  function,  Bq.  VI.6,  on  z  lo  the  correct  one,  the  coefficient  x^  can 
not  be  selected  so  ac  to  satisfy  the  integral  equation  In  the  whole  range  (z|  >  L. 
It  is  advantageous  to  select  a  coefficient  x^  which  satisfies  the  Integral  equa¬ 
tion  for  a  value  of  z  associated  with  a  relatively  large  value  of  ot(z)  and  herce 
with  a  large  contribution  to  the  far  field  potential,  viz.  for  z  *  L,  say  U  e: 

*1  *  “ 

/  r(Us-z')  •  (-72  +  ~  -75)  «p(iKzO  dz' 

L  2  *  2 


1!> 


(VI.1U) 


An  alternative  procedure,  which  gives  afire  nearly  equal  weight  to  the  whole  region 
of  z  where  the  integral  equation  apr  J.es,  is  based  on  the  fact  that  for  the  correct 
function  a(z),  J[a]  equals  -A [a],  fron  Sqs.  V.6  and  V.iQa.  Bence,  substituting  the 
ratio  x^/x^,  from  Eq.  VI.  13,  and  the  value  of  J[a]  from  Eq.  71.12,  in  the  expres¬ 
sion  for  A[o],  Eq.  VI. 9,  we  can  solve  for  x^ 


(VI. 15) 


Experience  with  numerical  calculation!)  will  indicate  which  procedure  is  preferable. 

So  refine  the  selection  of  the  trial  function  further,  we  can  introduce  addi¬ 
tional  unknown  coefficients  associated,  for  example,  with  naa-propagaulLg  inccet- 
prtssible  near  field  components  of  the  potentials .  The  trial  function  might  thus, 
for  exanple,  be  expressed  in  terms  of  three  coeff ic  Loots :  and  x^.  associated 

with  propagating  coaq>anents  of  the  potentials,  and  x^  with  an  incompressible,  near¬ 
field  component  decaying  rapidly  with  distance: 


(VI.6a) 


This  yields  three  simultaneous  equations  of  the  roro  of  Eq.  VI. 8.  Once  again,  ve 
will  find  that  tr»  se  equations  are  linear  ir  the  three  unknown  coeff icientf  end, 
of  course,  homogeneous.  Ve  can  therefore  construct  a  third  order  determinant 
similar  to  Eq.  VI. U.  The  constant  torn  and  the  linear  term  in  J|a]  ere  found  to 
cancel,  leaving  only  a  cubic  and  a  quadratic  tern  in  J(a].  The  detero inart  thus 
yields  a  single  root  J[a]: 

Jlaj  -  fcfOfeVfc)  +  a2<blb3’b?3)  +  83(blb2-4>  + 


+  2al02^b13b23'b3b12)  +  2  V3<b12b23*b2b13) 

+  3^b12b13_blb23^  *  ^blb2b3+2b12b13b23 

*  (bib23+b2b13+b3b12)rl  ( VI. 12a) 

We  then  solve  three  of  the  set  of  three  homogeneous  equations,  Eqs.  VI. 8,  for  two 
ratios  of  uaderterained  coefficients.  Finally,  we  solve  for  the  amplitude  of  the 
one  remaining  coefficient  by  satisfying  the  integral  equation  tt  s  «  I  ♦  «,  or  in 
the  manner  Indicated  in  Eq.  VI.15. 
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If  a(z)  is  expressed  in  terns  of  5  unknown  coefficients,  the  functionals 
A[a]  and  B[otJ  and  their  derivatives  take  es  the  following  fors: 


Alaj  =XVa  ’  sr  =  aa 


B[o]  = 


nSI 


2  )  b  x  x  ) 
/  mm' 


=  2b  x  +  2  )  bx 
n  n  /  .  ms 


(VI. 9a) 


The  set  of  K  homogeneous  linear  equations  for  thf  Jmovn  coefficients  corres¬ 
ponding  to  Eqo.  VI. 10  is  of  the  general  form 


2(a^bnJta])xn  ♦  2 


tanam*bm,JE0^X«  “  0 
do  an  n 


(Vl.lOa) 


When  the  coefficient  matrix  of  this  set  of  equations  is  set  equal  to  zero  It 
will  be  found  that  only  the  terms  containing  the  two  highest  powers  of  Jfa], 

W  1 

H  and  P-1,  do  not  cancel.  When  both  teres  ore  divided  by  (j[a))  ,  a  linear 

equation  in  J(cc]  is  obtained.  The  Nth  order  determinant  for  *T{cs]  has  a  single 
non-vanishing  root.  This  is  consistent  with  the  requirement  that  the  integral 
equation,  Eq.  IV. 10,  have  only  oca  solution. 

Experience  with  numerical  calculations  *111  show  whether  the  radiation  im¬ 
pedance  is  sensitive  to  the  selection  of  the  trial  function  a(z).  If  this 
should  be  the  case,  the  functional  dependence  of  the  near  field  potentials  on  z, 
particularly  of  the  non-propagating  incompressible  components,  can  be  studied 
more  closely  so  as  to  construct  a  more  sophisticated  trial  function  than 
Eqs.  VI. 6  or  6a.  Theoretical  insight  Into  this  functional  relation  can  bt>  gain¬ 
ed  from  the  fluid  mechanics  literature  dealing  with  accessions  to  inertia  of 
vibrating  solids.  Comparison  with  the  results  of  the  non- variational  solutions 
presented  in  Appendix  B  can  also  be  ussd  to  evolve  more  refined  expressions  of 
a(z). 

The  fact  that  J[a]  £r,  referring  to  Eq.  V.7,  Z^)  can  be  evaluated  from  a 
variational  principle,  without  previously  determining  the  amplitude  of  the 


A 


4 


1 


unknown  coefficients  of  the  trial  function  »(z),  has  already  been  related  to  the 
Larine -Schwinger  variational  principle  for  scattering  cross  sp'-tians  (Section  151). 
Another  parallel  which  say  be  more  familiar  to  svce  readers  is  found  in  the 
Rayleigh- Ritz  method  for  optimizing  the  natural  frequencies  obtained  from  Rayleigh's 
principle.  In  this  method  a  trial  function  is  assumed  for  the  dyrauaic  configura¬ 
tion  of  the  vibrating  system.  The  best  choice  of  the  coefficients  in  this  trial 
function  in  determined  by  giving  a  stationary  value  t.o  the  natural  frequency  ob¬ 
tained  from  Rayleigh's  principle.  If  B  unknown  coefficients  arc  used  in  expressing 
the  trial  function,  a  set  of  H  linear  homogeneous  equations  is  obtained  with  the 
coefficients  as  unknown  quantities.  By  setting  the  coefficient  matrix  of  this  set 
of  equations  equal  to  zero,  values  of  natural  frequencies  are  obtained,  without 
ever  having  to  compute  the  unknown  coefficients  themselves .  The  fundamental  natu¬ 
ral  frequency*  thus  obtained  is  equivalent  to  the  -unctiocal  Jfa).  To  ccaqmte  the 
ratio  of  the  undetermined  coefficients  at  that  frequency  one  substitutes  this 
value  of  the  fundamental  frequency  back  in  the  set  of  K  hcKogoneous  equations  and 
solves  for  H-l  ratios.  The  amplitude  of  the  Hth  unknown  coefficient  is  finally 
obtained  from  an  Inhomogeneous  equation  of  motion. 

To  conclude  our  study  of  the  solid  cylindrical  radiator,  we  now  turn  to  the 
evaluation  of  the  far  field  potentials. 

VII.  The  Fhr  Field  Potentials 

To  evaluate  the  potential  in  the  region  r  <  a,  z  >  L  us  substitute  the 
Green's  function,  G^,  Bq.  IV.2,  in  Eq.  IV,7a: 

»  °  J  (k  r) 

4i(r»z)  a  §jfa(s')dz?  J  {exp[lkz(z-z')1texp[lkz(n-t'-2L)]}dka.  (VII. l) 

J  *  t  1'  r 

L 

Fbr  a(z'),  we  substitute  s  trial  function  of  the  form  of  Rq.  VI. 6,  with  the  un¬ 
known  coefficients  expressed  in  terms  of  J[«],  as  described  in  Section  VI.  The 
k^,- Integral  associated  with  the  first  exponential  term  in  braces  in  Eq.  VH.  1  is 

*Xn  contrast  to  the  variational  principle  used  here,  which  yields  a  single 
solution  «(a] ,  the  Rayleigh-Rltz  techniques  yields  a  nuabvr  of  natural  frequencies 
ev’uai  to  the  mnier  of  unknown  coefficients  in  the  assumed  trial  function. 


is  given  in  Appendix  A,  E q.  A. 14.  Substituting  (-z'+2L)  in  place  of  s',  we  obtain 
the  kz-integral  associated  with  the  second  exponential  tern.  When  we  add  the  two 
integrals,  we  obtain 


dkz=2si[- 


(lklz-z 

ka 


Jo(knr) 


i  Z  ^  ^  fa 

[exp[i(k2-k2)^  jz-z'  |}+expte(k2  -k2)^(z+z'-2L)]}  (VH.2) 

where  (kQa)  is  the  nth  root  of  the  Bessel  function  The  integral  over  z'  in 
Eq.  VH.1,  aust  be  split  into  two  regions  of  integration:  (l)  from  L  to  z,  where 
Jz-z'l  is  taken  equal  to  (z-z')>  a»i  (2)  from  z  to  <*>,  there  |z-z'  |  equals  -(z-z'). 
Since  even  the  lowest  root,  k^a=3.83,  is  generally  larger  than  the  ka -value  of 
resonant  piezoelectric  or  aagnetastrictive  transducers,  the  terms  under  the  suazaa- 
tion  sign  decay  exponentially  with  increasing  Jz-z'l.  Because  the  source  distri¬ 
bution  a(z)  extends  to  infinity,  these  "near  field”  terms  contribute  to  the  far 
field.  By  using  energy  flow  considerations  it  was  shown  in  Section  VI  that  the 
desired  far  field  behavior  of  i. ,  Eq.  VI. 1,  requires  that  the  function  a(z)  embody 

i  i-2  1 

terns  of  order  |z|  and  hi^er.  She  dominant,  plane-wave  components  of  the  in¬ 
verse  transform  of  the  Green's  function,  Eq.  VH.2,  do  not  decay  with  increasing  z‘ . 
In  combination  with  the  far  field  term  of  cr(z'),  these  plane  wave  components  of 
the  Green's  function  therefore  give  rise  to  a  far  field  potential  whose  absolute 
value  varies  as  flz  1  2dz  =>  -  This  result  is  consistent  with  the  potential, 

Eq.  VI. 1,  used  in  deriving  the  functional  fora  of  a(z)  in  the  far  field.  The 
evaluation  of  the  inverse  transform  of  the  Green’s  function,  Eq.  A. 14,  is  thus 
.•oasis tent  with  the  energy  flow  analysis  in  &js.  VI. 2  to  4. 

The  far  field  in  the  region  r  >  a  is  obtained  by  substituting  the  appropriate 
Green's  function,  Eq.  XV.l,  in  Eq.  IV. 5s: 


CS>  3) 

"  r,  J  } 


Mm) 


exp[ik2(s-z/)]  u(z')  dk2  dz' 

r  / 


(ra.3) 


The  integration  over  z'  can  be  carried  out  immediately  by  making  use  of  the 
definition  of  the  Fourier  transform  of  the  velocity  distribution  u(z'): 


GP 

u(kz)  =J  u(z')  exp(-ik2z')  dz' 


(vn.4) 


The  expression  for  the  patentlcl  now  becoaes 
e  u(k  )H  (k  r) 

*o^)-fc  /  k&cvr  *. 


(vn.5) 


When  the  asynptotic,  large  argument  expression  for  the  Hankel  function, 

*"  (vn.6) 


Eo(V}  =  apUOyr-f}* 


is  substituted  in  the  integrand  in  Eq.  VII. 5,  and  using  spherical  coordinates  R 
and  $  in  lieu  of  the  cylinder  coordinates, 

2  =  R  cosS  and  r  =  R  sing, 

the  far  field  potential  be cooes 


a 

■  “rf-W*)  .  r 

(2*^R  Sin9)®  i. 


exp[iS(krsin9  +  k^cose)] 


snVkra) 


dk 


(vn.7) 


This  integral  vas  evaluated  by  laird  and  Cofcen*^  using  the  aethod  or  stationary 
phase.  The  result  thus  obtained  is 

•o<M>  ■  ('“•S» 

The  velocity  transfora  u(k  coss)  can  be  written  Ecre  explicitly  in  terras  of  the 
trial  function  a(z')  os 

u(k  cose)  »  Ja(  z‘)  cos  (tot'  eoss)  ds']  \VU.9) 


We  now  have  concluded  tbs  analysts  of  the  solid  cylindrical  radiator,  and 
proceed  with  the  open-ended  cylinder. 
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YIH.  The  Open-Ended  Free-Flooding  Cylindrical  Radiator  of  Vanishing 
Vail  Biictaiess 

We  consider  a  cylindrical  radiator  whose  wall  thickness  is  negligible  com¬ 
pared  to  both  its  radius  and  the  acoustic  wavelength.  For  such  a  source  configu¬ 
ration  we  can  construct  a  single  potential  for  the  region  r  <  a  extending  now 
from  -®  <  z  <  =>.  This  is  in  contrast  to  the  solid  radiator  where  we  had  to  define 
two  potentials  5  and  each  valid  in  a  seal-infinite  region.  The  outer  poten¬ 
tial  is  similar  to  the  outer  potential  derived  for  the  solid  cylinder  ,  2qs.  IV. 5a 
and  7a.  She  inner  potential  is  of  the  form 

CP 

^(rjz)  =  -2«a  f u(z')  G^r ,8,2-2*)  dz*  (VHI.1) 


where  the  Green's  function  G^  1»  <2»5  in  Eq.  IV. 2.  u-r^ring  this  with  Eq.  IV. fb 
we  see  that  the  potentials  defined,  respectively,  for  the  free -flooding  and 
solid  case  differ  as  to  the  range  of  s' over  whi-1  the  integration  is  performed  as 
well  as  to  their  Green's  functions.  As  in  the  case  of  Gq,  it  is  convenient  to 
define  separately  the  component  of  0^,  which  is  symmetrical  about  z'eO,  and  which 
alone  contributes  to  the  potential  when  the  velocity  distribution  is  similarly 
symmetrical: 


S^{v,a,2-2')  = 


(VHI.2) 


Assuming  a  constant  velocity  0  over  the  radiating  surface,  this  potential  is  again 
expressed  in  terns  of  an  unknown  velocity  distribution  0  3(2'): 

L  ® 

$A(r,2)  =  -  4xa  0[ J  gi(r,a,s-z')  dz'  +  Ja( z,)R1(r,o,z-z>)  dz']  (VHI.3) 

The  continuity  condition  at  the  boundary  r»a,  z  >  L,  is  a/^in  in  the  form  of 

Eq.  IV. 9.  The  corresponding  integral  equation  is  therefore  also,  formally  at  least, 

similar  to  the  integral  equation  derived  for  tne  solid  cylinder,  Sq.  IV. 10.  Hov, 
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however,  the  kernel  T(z-z')  and  the  non-hooogeneous  tern  2(2}  are 
r(z-z')  »  ^(s.ajz-z')  +  gQ(e,a,z-z') 

L 

H(z)  =  j  [g^a.a^-z')  +  8o(aJa,2-2/)]  dz' 


(vm.te) 

(VIXI.Ub) 


The  expression  for  r(z-z')  can  he  siaplified  by  using  the  Wrooskian  relation  for 

1  f* oM  vv) 


H  and  J  : 
o  o 


.  r  HIM)  J  1X8} 
=  ^ 1 


cos  k  2  co6  k  z  dk 
z  z  z 


i  f  cos  kza 
*3a2J  k?J,(k_a 


cos  k.z  cos  kzz* 


0  m'¥ra> 


dk 


(vm.5) 


The  radiation  iapedance  of  the  free- flooding  shall  differs  froa  that  of  the  solid 
cylinder  in  that  the  pressure  on  both  the  outer  and  inner  surface  contribute  to  it 

L 

(VHI.6) 

v  * 

\ 0 


u 

2  .  .  *52_£2£  J  [ya.z)  -  ^(a.z)]  dz 


The  tvc  cospcoents  of  this  iapedance  stated  in  Eq.  IH.1  cen  again  be  separated. 
The  iapedance  associated  with  Robey’s  ne these tical  nodel  is: 


ZT  =  -(bsa)2  top  j  J  r(z-z') 


dz'  dz 


0  0 
L 

p  r 

■  -(!<aa)  topj  H(z)  dz 

The  correction  tens  resulting  free  fluid  flow  across  the  cylindrical  boundary 
r=a,  |z  |  >  h  is 

L  • 


(vra.7) 


Za  =  -(bsa)2  top  J" [J a(z')  r(z-z')  dz']  dz 


(VIII. 6) 


0  L 
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Since  T(z-a/)  is  syintetrical  in  2  ami  z‘ ,  vs  can  invert  the  order  of  integration 
and,  using  the  definition  of  H(z)  in  Bq.  VlH.Ub,  write  the  uninova  coeponent  of 
the  radiation  isspedance  as  follows: 


ZQ  =  -(bita)2  top J  w(z)H(z)  dz 


fVTTT  a\ 
V  —  -  '  » 


The  construction  of  the  functional  and  the  proof  that  it  is  stationary  for  the 
correct  fora  of  the  unknown  function  a(z)  parallels  foraally  the  proof  given  in 
Eqs.  V.8  to  16  for  the  solid  cylinder.  The  relation  between  the  unknown  impedance 
and  the  functional  j[a],  Eq.  V.7,  is  also  applicable.  The  variational  solution 
of  the  free-flooding  thin-walled  shell  is  therefore  foraally  identical  with  thet  of 
the  solid  cylinder  provided  the  definitions  of  r(z-z')  and  H(z)  given  in  Eqs.  VUI.ha 
and  b  are  used,  instead  of  the  corresponding  definitions,  Eq.  IV. 11  and  8,  re¬ 
spectively,  which  apply  to  the  solid  cylinder  case .  Ve  will  see  in  the  next  sec¬ 
tion  that  this  parallel  does  not  hold  when  ve  assume  a  realistic  free-flooding 
transducer  or  "squirter"  whose  wall  thickness  is  not  negligible. 

The  expression  for  the  far  field  potential  9  is  still  given  by  Eqs.  VH.8  and 
9.  The  expression  for  the  potential  is  somewhat  different,  because  of  the  con¬ 
tribution  of  the  region  |  z'|  <  L  which  is  absent  in  the  case  of  the  solid  cylinder: 

.  r  C  3  (kr) 

Si(r'z)  =  fe  J  u< z>)  J  e3CP^z(*-z')J  dz'  (VIII.  10) 

The  integral  over  dk^  is  given  in  Appendix  A,  Eq.  A. It.  The  consents  made  in 
connection  with  Eqs.  VII. 1  and  ?  apply. 


EC.  The  Free-t -oodl 


ylindrlcal  Transducer  or  "Squirter" 


When  we  drop  the  assumption  of  a  vanishing  wall  thickness,  we  must  formulate 
the  analysis  in  terns  of  three  coaxial  cylindrical  boundaries  and  their  respective 
radial  velocities  (Pig.  3). 

(l)  ^^(l  -  £);  radial  ve;"  ty  u^z) 


(2)  r=a 


;  radial  velocity  u(z) 


(3)  r=a  =a(l  +  radial  velocity  u  (z) 


(K.  1) 


■ 


Each  of  these  cylindrical  boundaries  is  of  infinite  extent  In  the  a -direction. 

The  required  modifications  in  the  expressions  for  the  potentials,  Eqs.  IV. 5a 
and  vm.l,  are  self-evident  and  Involve  merely  labeling  a  and  u  with  the  appro¬ 
priate  subscripts  i  or  o.  There  is  an  equally  obvious  change  in  Eq.  VIU.6  for 
the  radiation  tape  dance,  where  the  potentials  and  #  must  now  be  multiplied, 
respectively,  by  the  ratios  a^/a  and  a^/a.  A  noa-trivicl  change  must  be  intro- 
duced  in  the  statement  of  the  continuity  condition,  which  now  no  longer  takes  the 
fora  4j(a,s)=So(a,s).  Bather,  we  as suae  an  Incompressible  potential  in  the  annu¬ 
lar  i\  gioo  ai  <  r  <  aQ.  With  this  assumption  the  difference  between  the  poten¬ 
tials  5.  and  $  must  be  aatched  to  the  inertia  force  exerted  by  the  fluid  located 
in  this  annular  region. 

(l  +  *0(VX*  *  a  -21C>(z)U,  for  jz|  >  L  (K.2) 

The  assuaqition  of  an  incompressible  potential  implies  that  the  ratio  2h/x  is 
snail.  This  condition  is  generally  satisfied.* 

We  will  now  derive  relations  between  the  radial  velocities  on  the  three,  cylin¬ 
drical  surfaces  defined  above.  As  a  result  of  the  assumption  of  an  incompressible 
potential  in  the  annulor  region  ai  <  r  <  eQ  the  velocities  in  the  two  semi- 
infinite  regions  prolonging  the  transducer  can  be  derived  from  the  requirement 
that  inflow  must  balance  outflow  across  the  cylindrical  boundaries  r*«,  o^,  and  aQ 

u^  »  ua/aj^ 

°  u/(l-h/a) 

U  a  uo/a 
o  o 

a  u/(l  t  h/a)J  (DC. 3) 

*For  the  !TRL  oaguetcctr tetire  transducer  ring  2h/\  is  approximately  0.07*  For 
transducers  where  this  assumption  is  not  valid,  a  compressible  potential,  4Q,  must 
be  constructed  for  the  region  <  r  <  a^.  The  three  potentials  must  oetisfy  two 
continuity  conditions,  viz:  *i(°i),!®a(8i)  aai  Sa^8o^=*o^0o^’  ^nsteaii  of  one>  tvo 
unknown  radial  velocity  distribution  over  the  boundaries,  r»a^  and  r*ao  coot  be  de¬ 
termined  from  the  two  simultaneous  Integral  equations  arising  from  the  two  poten¬ 
tial  continuity  requirements.  These  equations  have  Been  constructed,  but  It  has 
not  yet  been  verified  whether  a  variational  principle  can  be  applied  to  their 
solution. 


For  Izj  <  L,  i.e.,  in  the  region  of  the  transducer,  different  relations  must  be 
used.  The  fom  of  the  equations  relating  these  three- velocities  depends  upon 
whether  we  are  dealing  with  piezoelectric  or  aagnatoatrictive  elements.  In  the 
case  of  piezoelectric  transducers,  the  vcltege  applied  to  the  electrodes  located 
on  the  inner  and  outer  surfaces  of  the  ceramic  ring  produces  a  radial  strain  *_. 
The  corresponding  eircu  iferential  strain  «  of  the  mean  surface  of  the  element  is 
t>«»  result  of  Poisson  coupling: 


e_  =  -  w  e 
<p  r 


where  v  is  Poisson's  ratio.  This  circumferential  strain  is  related  to  the  radial 
velocity  u  of  the  mean  surface  as  follows: 


€ai  =  a/(-i®a) 


(DC.5) 


Hoting  that  e_  is  of  opposite  sign  than  c  and  hence  than  the  displacement 
-1  ^  ^ 

u(-lxo)  ,  we  find  that  the  velocitv  of  the  outer  and  inner  surface  are  recpective- 
ly  reduced  and  increased  by  the  radio i  strain 

=  u  -  hcr 


u  =  u  +  he_ 
o  r 


(IX.6) 


where  a  contraction  corresponds  to  negative  tr-  Combining  these  equations  we 
finally  have 

ui 3  u<x  + 


uo  ”  u(l  * 


(K.7) 


In  the  case  of  a  ring-shaped  laBgnetostrictive  transducer  the  current  in  the 
solenoid  produces  a  circumferential  strain  c  .  In  this  case  it  is  the  radial 
strain  that  results  from  Poisson  coupling: 


=-v  e 
r  :p 


v  u(-ima)” 


(IX. 8) 


Combining  these  equations  we  now  have 

u.  =  +  — ) 

u  8 


u0  3 


(XX.9) 


it 


For  the  sake  of  brevity,  we  define  the  following  coefficients,  which  tend  to  unity 
for  snail  values  of  h/a: 


»j  -  <i »  ;)  &  t  f) 


for  piezoelectric  transducers 

for  aagnetostrictive  transducers 


(DC.IO) 


If  ve  now  substitute  these  velocities  in  the  expressions  for  the  potentials  used 
in  the  boundary  condition,  Fq.  IX. 2,  we  obtain  a  acre  complicated  integral  equa¬ 
tion  than  for  the  two  earlier  configurations : 


CD 


(1-  |)81(*1»a.i,2-z,)3  <&' 


-I 


C80(l+  |)g0(a0»a0>s-2')  +  BjU-  |)s1(a1,a1,z-s/)]  dz*  forjz  j  >  L 


(IX.U) 


As  will  be  seen  shortly,  the  presence  of  the  linear  fc-r  m  a{z)  outside  the 
integral  sign,  which  aaXcs  this  into  a  Fredholm  1st  ."l  equation  of  the  second 
kind,  does  not  interfere  with  tie  applications  ef  V.  varl  tional  principle.  The 
presence  of  the  coefficients  and  Bq  in  the  aon-hosogeneou s  term  of  the  integral 
equoticas  does  unfortunately  sake  the  application  of  the  variational  principle  in- 
practical  because  the  stationary  potential  J(r)  which  (an  be  constructed,  is  no 
longer  proportional  to  ZQ,  as  stated  in  Eq.  V.J.  To  sake  this  single  relation 
applicable,  we  must  assure 


8*0  SI 
I  o 


(DC. 12) 


The  error  in  this  procedure  is  seen  froo  Eq.  IX *10  to  be 


e 


i  «*•>  - 


(DC.  13a) 


for  eognetostrictive  transducers.  For  piezoelectric  transducers  the  error  is  larger 


(DC. 13b) 


She  staple  variational  technique  is  therefore  better  suited  to  nagastostrictive 
than  to  piezoelectric  transducers.  For  the  nagnetostrictive  NRL  transducer 
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ring  (2a=5-7/8  la.,  2h=£  la.),  the  error  computed  frees  Eq.  EC.  13a  la  approximately 
6  percent.  It  may  seen  inconsistent  to  Introduce  ‘Assumption  DC.  12  and  not  to  drop 
the  linear  term  in  a(z)  from  Eq.  DC, 11,  and  the  ratio  h/a  from  terms  of  the  form 
(l  t  h/a).  Further  work  is  necessary  to  determine  whether  retention  of  these  terms 
increases  accuracy.  Until  this  is  done,  we  shall  retain  these  terms,  because  they 
do  not  complicate  the  variational  technique.  In  addition  to  introducing  the  assump¬ 
tion  stated  in  Eq.  DC.U,  we  give  a  new  definition  of  the  function  r(s-z')  and  of 
H(z) 

r(z-2')  =  (l  +  |)g(j(ao,oo,z-z/)  +  (l  -  -Jg^a^a^z-z')  (BC.llw) 

L 

H(z)  =  J  t(l  +  £)g3(aa,a0,z-z')  +  (l  -  })et(«1,a1,B*'))  <**'  (DC.l4b) 

Furthermore,  we  make  the  integral  equation,  Eq.  DC.U,  formally  into  a  Fredholm 
equation  of  the  first  kind  by  the  artifice  of  adding  a  Dirac  delta  function 
5(z-z')  to  the  kernel: 

CD 

J  [r(z-z')  +  6(z-z')]  a(z')  dz'  =  -  H(z),  for  |z|  >  L  (Dt.15) 

h 

The  functional  J[a]  is  still  of  the  seme  form  as  in  the  two  earlier  analyses, 

Eq.  V.9,  and  the  definition  of  the  functional  A[a],  also  regains  formally  the  same, 
Eq.  V.lOa.  The  functional  B[a]  is,  however,  different. 

®  cc 

B[a]  *  J  j  a(z)[r(z-z')  +  Jlj  «(z-z')]  a(z')  dz  dz'  (K.l6) 

L  L 

Ve  will  now  show  that  the  stationary  character  of  Jfcr]  can  be  established  as  before. 
Setting  the  increment  of  the  functional  equal  to  0,  we  have 
® 

6j[a]  a  o  a  j  a(z)  6a(z)  dz  - 

L 


-  4^  f  f  [r(z-z')+  ~6(z-x')Ka(z)ta(z')««(z')6a(z)ldz  dz'  (IX.  17) 

B2[a]  l  J  ** 
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The  order  of  integration  of  the  0(2)40(2')  tern  in  the  double  integral  can  oe 
inverted. 


CO  CD  CO  O 


// 


^  J  sff  [r(a-z')  +  Ha-z*)]  a(z')  4a(z)  dz  dz' 

LL  L  L 


(DC.18) 


From  here  on,  the  proof  parallels  exactly  the  steps  fron  Eqs.  V.13  to  16  ana  vlil 
therefore  not  be  repeated.  The  radiation  impedance  correction  factor  ZQ,  is  once 
again  formally  given  by  Eq.  Vm.S,  with  r(z-z')  defined  in  Eq.  E.lha.  Thus  by 
setting  the  coefficients  0  equal  to  unity,  ZQ  can  still  be  expressed  in  terns  of 
the  functional  J[a]  as  in  Eq.  V.7.  The  component  of  the  radiation  Impedance  asso¬ 
ciated  vlth  Robey’s  mathematical  model  is 

L  1 

V  -i»o(hw)2|  J [Bo(l+  |)go(ao,ao,z-z')+81(l.  fjg^a^z-z')]  da'  dz  (BC.19) 

Even  if  the  coefficients  0  had  not  been  set  equal  to  unity,  a  stationary  po¬ 
tential  J[a]  could  nave  been  constructed  vlth 

A[a]  *  J  a(z){  j '*  (B0(l*  f}go(a0,a0,z-*')+3i(l-  jj)^®^.*-*')]  dz'Jdz  (C.20) 


instead  of  the  expression  in  Eq.  V.lOa.  The  usefulness  of  the  variational  method 
is  hovever  impaired,  because  the  radiation  impedance  component  Z^  does  not  change 
in  the  same  manner  as  A[o]:  Za  13  given,  as  before,  by  Kq.  Vm.8  with  F(z-z')  as 
defined  in  Eq.  jDC.lUa.  It  therefore  does  not  involve  the  coefficients  0,  whether 
ve  3et  the  coefficients  equal  to  unity  or  not.  ZQ  is  therefore  not  proportional 
to  A[a],Sq.  IX. 20,  and  Eq.  7.7  relating  Za  and  J[o]  does  not  apply.  Thus  even 
though  the  vorlatloni-1  method  can  be  used  for  "squirters"  whose  walls  are  too  thick 
to  permit  setting  the  coefficients  S  equal  to  unity,  the  unknown  coefficients  in 
the  trial  function  a(z)  carat  bo  solved  for  before  computing  ZQ.  Whether  such  a 
procedure  is  competitive  with  the  non -variational  solutions  presentee  in  Appendix  B 
for  thick-walled  "squirters"  can  best  be  verified  empirically  after  numerical  cal¬ 
culations  have  been  performed. 

We  shall  now  extend  the  variational  technique  to  arbitrary  non-axisyaaetric 
valocity  distributions  of  the  radiating  surface. 


flinders  Vibrating  in  Longitudinal  and  la  Hon-Axlsyaaetric  Hades 

The  analysis  of  the  "squirter"  and  of  the  solid  cylinder  can  be  directly  adapted 
to  the  case  of  nonuniform  axisymaetric  velocity  distributions  over  the  radiating 
surface,  by  replacing  the  constant  velocity  U  in  the  region  |z |  <  L  vith  a  z- 
dependent  velocity  u(z).  If  u(s)  is  not  syaaetrical  about  z=0,  the  most  convenient 
approach  is  to  consider  the  velocity  distribution  as  the  sua  of  a  syaaetrical  dis¬ 
tribution  ug  and  of  en  antisymmetric  distribution  ug.  As  we  are  dealing  with  a 
linear  problem,  ve  can  add  the  corresponding  potentials.  We  first  confute  the 
potential  associated  with  u  as  in  the  preceding  analysis,  setting  at  (z)=a  (-z), 
and  using  the  corresponding  Green's  functions  g0>  Eq.  IV.6,  and  g^,  Eq.  VIII. 2 
(the  latter  in  the  case  of  the  open-ended  cylinder).  To  this  we  add  the  potential 
resulting  from  the  velocity  distribution  u  for  whi-h  a  (z)  »  -a  (-z).  The  suit- 
able  partial  Green's  functions  are  obtained  by  modify 'eg  tie  expressions  for  gQ 
and  g1  given,  respectively,  in  Eqs.  IV. 6  and  VIH.2,  sin  A^z  sin  k^z*  being  sub¬ 
stituted  for  the  product  of  cosines.  We  must  thus  solve  two  uncoupled  Integral 
equations  for  the  two  unknown  velocity  distributions,  a  and  a  .*  Unless  we  pro- 
ceed  in  thi3  fa p hi on,  the  solid  cylinder  with  arbitrary  velocity  distribution  u(z) 
gives  rise  to  three  different  potentials  *i+,  which  in  turn  result  in  two 
distinct  boundary  conditions  corresponding ,  respectively,  to  the  regions  z  <  L  and 
z  >  L.  The  two  resulting  integral  equations  will  thus  be  coupled,  each  involving 
both  unknown  velocity  distributions,  a(z  <  0)and  a(z  >  0). 

In  the  case  of  a  piston  or  ring  vibrating  in  phase  on  a  finite  cylindrical 
baffle  or  array,  the  velocity  distribution  of  the  active  element  is  of  course  con¬ 
stant  and  hence  syszaetrical  over  the  midp  lane  (zi=0)of  the  element,  but  unless  this 
element  is  centrally  located  with  respect  to  the  baffle  or  array,  the  velocity  dis¬ 
tribution  a(z)  will  not  be  symmetrical.  In  this  respect,  the  present  mathematical 
model  differs  from  Robey's  model,  in  which  the  potential  and  the  velocity  distri¬ 
bution  ere  always  symmetrical  about  the  plane  of  syasetry  of  the  active  element. 

A  configuration  of  practical  interest  is  that  of  a  wolid  cylinder  whose  end 
caps  reciprocate  in  the  axial  direction.  This  situation  arises  as  a  result  of 
Poisson  coupling  vith  predoalnantly  radial,  axisyasitrlc  modes.  End  cap  motion 
can  contribute  the  oajor  portion  of  the  sound  field  in  the  case  of  the  so-called 


*This  procedure  will  be  illustrated  in  the  report  dealing  vith  an  array  of  ring 
transducers  (see  footnote  on  p.  l). 
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accordion  aodes,  which  are  predoaimntly  longitudinal.  To  account  for  an  axlsyn- 
netric  velocity  distribution  v(r')  over  the  end  caps,  the  following  integral  is 
added  to  the  surface  integrals  in  Kq.  IV. 5b: 


(r,z) 


^ajv(r')  Gt(r,r' 


,2  +  l)  dr' 


(X.l) 


where  v  has  been  taken  positive  in  the  positive  z-directioo.  If  the  cylindrical 
surface  of  the  radiator  is  aotionless,  the  potential  $q,  in  the  region  r  >  a  is 
associated  entirely  with  the  velocity  distribution  a(z)  across  the  two  surfaces 
(r=a.  z$L). 

B»e  variaticoal  analysis  ;an  be  extended  further,  to  include  an  arbitrary  non- 
axisynnetric  velocity  distribution 


u(v,z) 


VB(z) 


for  J  z  l  <  L 


(X.2) 


is  a  nodal  velocity  amplitude,  tbe  naxinua  value  of  the  function  fa(z)  being 
unity.  To  each  Fourier  component  Da  of  the  velocity,  corresponds  a  partial  po¬ 
tential  »a(r,z),  the  total  potential  being  of  the  fora 


*a(r,z)  cxp(imp) 


(X.3) 


The  partial  potentials  ore  obfcf.ir-  ;  trass  Eqo.  IV. 5  and  Eq.  VUI.i,  by  substitu¬ 
ting  in  place  of  the  axisywntrlc  Green’s  functions  given  In  Eqs.  IV.  1  and  2  the 
following. 


1  f 

z-z')= --~-cxp[lB(ts^<p')]  f  -  ^ygp[ikz(z-zr))dhs>  f or  rbe  (x.h) 

Jf  S  T"  Z  * 

The  near  field  value  of  this  integral  has  been  evaluated  by  Greenspan  and  iherwan.® 
Its  asycptotic  far  field  valus-  is  gives  by  Laird  and  Cohen: 

■  if?S;  •  -  — 8  (*-5> 
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The  corresponding  Green's  function  in  the  cylindrical  region  Is  derived  In  Appendix 
A,  £$•  A.6: 

,  f  Jjkr} 

<5<1f,(r^»2-s#,<?-<p')=  ~“~exptia(f-V)]  I  r-ynft-y  for  Ks  {X. 

**  k/a  V/r*'  z  * 

•CD 

The  Integral  Is  evaluated  in  Eq.  A.l3.  Each  potential  4^  is  the  sua  of  two  cesspc- 
nents:  a  potential  i _ associated  with  the  knows  modal  velocity  distribution 

tSP 

Of  (i)  over  the  radiating  surface,  and  a  conpocent  $ _  associated  with  the  unknown 

g  a  OPC 

modal  velocity  distribution  0  a  (z)  in  the  two  regions  iz{  >  L-  The  former  sob- 

2*  d  * 

ponent  ^  is  of  course  tins  component  computed  fret  Robey’s  aetheoaticel  aodel  of 
tive  cylinder  prolonged  by  two  seal-infinite  cylindrical  baffles.  She  modal  imped¬ 
ance  associated  with  the  nth  node  of  the  radiator  c®  be  expressed,  as  in  Eq. 111.1, 
as  the  sua  of  hooey's  isgedancc,  associated  vi*-b  #  ,  and  of  a  correction  term 

associated  with  i  : 


4J 


i  ’nr1 

-l 


U 

z«v =  -  I  os 


This  isqjedance  can  hs  used  to  ecspute  the  generalized  force  associated  with  radia¬ 
tion  loading  of  the  Eth  elastic  safe  of  the  cylinder,  and  heuee  the  modal  imped¬ 
ances  snd  natural  frequencies  of  the  subsss  god  cylinder.  Modal  radiation  imped¬ 
ances  cen  also  be  combined  to  compute  the  self-radiation  impedance  of  rigid 
pistons  in  finite  cylindrical  baffle*.*  Because  of  the  similarity  in  the  fora  oS 
the  Green's  functions  of  ths  axisymaetrie  case  analyzed  in  detail  in  this  report 
and  of  the  non-azlsytssstrlc  radiator  configurations,  it  is  obvious  that  the  inte¬ 
gral  equations  which  aQ(z)  wust  satisfy  are  of  the  same  form  as  the  Integral 
equations  -ditch  define  a(z)  in  the  axisynsetric  ease.  Functionals  *TEIaBl  station¬ 
ary  with  respect  to  the  correct  function  aQ(c)  can  be  constructed  and  are  found  to 
be  of  the  same  form  as  ti  functional  J[cj]  constructed  earlier  for  the  exlsycset- 
ric  radiator.  55s  proof  that  the  iepedance  Z^  is  proportional  to  Jn(dn]  for  the 


*The  component  Z^  of  this  Lcpe  dance  is  computed  by  Greens poa  and  Sherman- 
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V.  *  <V  ?•  . 


cornet  fora  of  £_(;■}  para  Ole  Is  the  proof  in  Section  7.  3he  variational  technique 
*» 

illustrated  in  Section  VI  sen  therefore  be  used  to  evaluate  2^,,,  and  need  not  te 

jgX' 

repeated  here. 

Is  the  case  of  tfcs  solid  cylinder,  ncarigid  vibrating  end  caps  can  be  accounted 
for  by  adding  to  the  expression  for  the  potentials  ft  surface  integral  over  the 
eca  caps: 

*  Jf  a 

45^+  =  +  2  j  |  v(r',9')  G^C  r,r',  <^',2  +  L)  r"  dr'  d©‘‘  (X.8) 

IS=li  0  (> 


Appendix  A* 


BSBIVAIIOB  AHB  EVA1CAKCS  OF  TBS  QHStJt'S  FUHCTXCB  G. 

FOR  THE  CYLIKDRIGAIi  RSGIOS  r  <  a  1 


1.  Construction  of  the  Green's  Function 

This  derivation  parallels  the  construction  of  the  Green's  function  G  for  the 

6  0 

region  r  >  a  given  by  Robey  for  Eeuaenn  boundary  conditions  (d*/dn  known,  <3G/de 

«  Q 

made  to  vanish  on  boundary)  and  by  Papas'10  for  Dirlchlet  boundary  conditions  ($ 
known,  G  node  to  vanish  on  boundary).  The  Green's  function  associated  with  the 
ath  Fourier  cceponent  of  the  velocity  distribution  in  <p  can  be  expressed  in  terns 
of  an  inverse  Fourier  transform  in  (s-z'): 


f  Gla(r»r/,^z)exptik3(2-a/))  dkz 


The  Green's  function  satisfies,  by  definition,  the  non-hoeogeneous  Helmholtz 
«w>ation  expressed  in  cylindrical  coordinates.  Consequently,  the  transform  of  the 
Green's  function  satisfies  the  following  equation: 

*f  I?  {r  +  **  *  Gm<r'r''kt>  "  '  (A’2) 

r 

2hs  step.'  leading  from  the  nac-hoeogenesnis  Helmholtz  equation  to  Eq.  A. 2  are  pre¬ 
sented  in  detail  in  reference  U.  Except  when  r=r/,  Eq.  A.£  is  of  the  form  of 
Bessel's  equation.  A  suitable  solution  to  this  equation  must  be  regular  when 
r  <  vanishes  (r  <  and  r  w  are,  respectively,  the  scalier  and  the  larger  of  the 
quantities  r  and  rf).  The  solution  of  Eq.  A. 2  oust  therefore  contain  only  Bessel 
functions  of  argument  proportional  to  r  <  ;  Heuman  or  Henkel  functions  can  only 
hsve  Grouse nt 3  proportional  to  r  >  or  e.  A  combination  of  cylinder  functions 
which  satisfies  these  conditions,  and  whose  radial  derivative  dG/dr*  vanishes  on 


This  material  is  included  here  because  it  does  not  appear  to  be  available  in 
the  literature.  Appendix  A  is  condensed  from  a  Harvard  Acoustics  Research  labora¬ 
tory  Memorandum.1  r  Analytical  details  and  proofs  which  had  to  be  omitted  to  keep 
the  length  of  this  report  within  reason,  can  he  found  in  reference  1?. 


:v: 


■v.  ' 


the  cylindrical  boundary  r  >  =  r\  is  of  the  fora: 

®ia(r*r*  *  Va*-V  <)tH'(V)jB»krr  >}  *  Bm(V  >)ja^ra)]  'A 

The  coefficients  A^  are  determined  from  the  equation  defining  the  discontinuity 
of  the  first  derivative  of  the  Green's  function: 


r  5  yr/,  k8)  | 


After  sees  transformations  the  coefficient  A  is  found  to  be 

a 

*■ '  <1-5) 

When  we  substitute  this  coefficient  in  Eq.  A. 3,  set  r*  »  r  >  =  a  and  r  =  r  <,  and 
use  the  Wronskian  relation  between  H^Ck^a)  and  J  {k  a),  we  finally  obtain  the 
following  expression  for  the  Green's  function 

Gia(r,a-<S>V>3-z/)  “  HElidpSeJI  J  ^  «q>[lkz(g-n#)]  dkg  (A. 6) 


For  the  axisymetric  case,  we  have  vfaich  the  Green’s  function 

given  in  Eq.  IV. 2.  We  will  follow  the  notation  used  in  the  body  of  the  report, 
whereby  the  subscript  m  is  emitted  when  s=0,  i.e. ,  in  what  follows ,  G^,  I,  Rq  and 
kn  indicate,  respectively,  0iQ,  IQ,  and  k^. 


2.  Evaluation  of  ".he  Inverse  Fourier  Transform 

We  will  not  evaluate  the  infinite  integral  in  2q.  A. 6.  For  the  purpose  of 
analysis  the  wave  number  is  assumed  to  have  a  snail  imaginary  component 

ks  =  5  +  in  (A. 

The  integration  in  Eq.  A. 6  will  be  performed  in  the  complex  plane  along  a  closed 
counter clockwise  contour  Including  the  real  arts  and  a  half  circle  of  Infinite 
radius  (Fig.  U).  By  the  residue  theorem,1^  the  value  of  the  contour  integral  is 
2a i  times  the  sum  of  the  n  residues  at  the  poles  k_  of  the  integrand. 


^  . 


<0 

J  IBdkz  =  2«i  ^  ,  for  2-2'  >  0,  and  in  upper  half -plane 
-®  n 

V  • 

=”2Si  A_Rb21*  *  f°r  Z”z/  <  °  ’  and  ^nn  10  lover  half-plane  (A. 6 

n 

■v&ere  1  stands  for  the  Integrand  in  Eq.  A. 6.  This  contour  integral  equals  the 
integral  along  the  real  axis  if  the  contribution  of  the  half-circle  vanishes.  To 
achieve  this  condition  the  integrand  oust  vanish  as  kz  tends  to  infinity,  i.e., 
jexp(ik,(z-2/)}|  in  the  integrand  oust  decrease  exponentially  with  increasing  k^. 

Hence 

1}  >  0,  for  z-z*  >  0  (integration  in  upper  half -plane) 

*j  <  0,  for  z-z'  <  0  (integration  in  lower  half -plane)  (A. 9 

Like  kz,  k  can  be  assured  to  be  complex.  Its  infinite 3 iael  imaginary  component 
can  be  associated  with  viscous  losses  in  the  acoustic  medium,  if  a  physical  in¬ 
terpretation  is  desired.  The  complex  quantity  +k  lies  Just  above  the  real  axis, 
and  -k  just  below  it.  The  contours  of  Integration  are  then  as  shows  in  Fig.  U. 

He  will  first  evaluate  the  axlsymetric  Green's  function.  The  integrand  I 
has  poles  at  kr=0,  i.e.,  at  k2a  t  k.  Taking  the  asymptotic  expression  of  the 
Bessel  functions  JQ  and  as  their  argument  tends  to  zero,  we  find  that  the 
Integrand  tends  to 


2  e*p[lkz(k  -z')! 

X(v  ) - - - £—5 -  ,  as  k  -  t  k 

2  a2(kz+k)(kz-k) 

The  two  sinpie  poles  at  k^-  t  k  give  rise  to  the  following  residues 


,  for  2-z '  >  0 


(A. 10) 


,  for  z-z'  <  0 


(A. ll) 


Other  poles,  all  of  then  simple,  occur  at  the  roots  (kna)  of  ^(k^a).  The  corres¬ 
ponding  residues  are 


±2  exp[4i(k2-kM(2-s')]  Jn(k  r) 


*V-k2)3 


(A. 12) 


Bt -cause 


R^Xz-s')  =  -Rn.(z'-z)  , 


(A. 1’) 


the  integral  can  be  stated  as  follows,  without  regard  for  the  relative  Magnitude 
of  z  and  z/ : 


2*irexn(lklz-z'l)  .  ^^[i^-h^lz-zl}  J0(y) 

a  1  ha  +  7  .  ,„2  k2^.  .  0U(k_a)-J„( 


(h 


(A.1U) 


For  all  roots  k  a  of  j'  which  exceed  ka,  the  terms  under  the  suaaation  sign  decay 
n  u  l  . 

exponentially  with  increasing  |z-z,|.  For  higher  order  roots,  the  terms  under 
the  suswtion  sign  are  proportional  to 

ggLg&bA),  for  n  large  (A.l5) 

n(ra)* 

The  series  expression  in  5q.  A.lU  is  thus  seen  to  be  convergent  except  for  z=z*, 
which  fulfills  the  require  be  nt  of  a  Green's  function. 

The  convergence  of  the  Green's  function  as  j z  |  tends  to  infinity  is  not 
spherical,  but  relies  on  the  small  imaginary  component  of  the  wave  number  k.  The 
reason  is  that  we  have  constructed  a  Green’s  function  suitable  for  cylindrical 
region,  viz.  a  circular  pipe,  where  only  viscous  losses,  but  no  spreading  losses 
occur.  The  potential  ?i  does,  however,  vary  a3  I  zl"1  exp(ikz)  for  large  |z|, 
because  the  radial  velocity  a(z)  gives  rise  to  a  net  outflow  of  acoustic  energy 
from  the  region  r  <  a  (see  Eqs.  VI. 1  to  6,  and  consents  following  Eq.  VH.2). 

We  now  turn  to  the  evaluation  of  the  non-axlsymbetric  Green's  function.  At 
kz=  ±  k,  for  a  >  0  the  integrand  tends  to 


I  _  -fa/?!..  exp[ik  (z-z')J  as  k  -  ±  k 

OB  Z  Z 


(A. 16) 


There  are  therefore  no  poles  at  k2°  t  k,  only  for  n=0.  For  njfo,  all  the  residues 
are  associated  with  the  roots  of  J*: 

Q 


w 


texp[tit^-k^.)‘(z-z'))  Jm(hmr) 

a2(k2-k2  )* 
cn 


for  njfo 


(A.  17) 
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Once  again  the  integrand  can  he  expressed  vtthout  regard  for  the  sign  of  (z-z')t 

for  m/0  (A.  18) 

Ibe  higher  order  terns  are  again  found  to  he  proportional  to  the  expression  in 
Eq.  A. 15,  and  thus  to  conform  to  Green's  function  requirement  hy  convertring  for 
zjlz'.  A  proof  was  civen  in  ref.  17  of  the  fact  that  even  though  the  function 
(k  -k„)*  has  a  branch  cut  in  the  region  |k2 |  <  k  of  the  real  axis,  the  integrand 
in  Eq.  A. 6  does  not  have  branch  points  at  any  of  its  poles. 


37 


Appendix  B 

HOH-VAHIAHOHAIi  T35CHHKWES  FOB  SOLVUG  THE  "SQUEffiER"  DKEGRAL 
EQBATIOH,  EQ.  EC.  11 

In  Section  EC  it  was  shown  that  the  variational  technique  developed  for  the 
solid  cylinder  is  applicable  to  the  free-flooding  cylinder  only  when  the  coeffi¬ 
cients  snd  &o,  Eq.  DC. 10,  can  be  set  equal  to  unity,  i.e.,  when  the  ratio  of 
wall  thickness  to  radius  2h/a  is  stall.  It  was  also  pointed  out  that  when  the 
ratio  of  wall  thickness  to  acoustic  wavelength  2h/X  is  not  snail  enough  to  make 
the  cspapressibility  of  the  fluid  annulus  in  the  region  <  r  <  n^  negligible,  a 
cooplicated  analysis  involving  two  coupled  s iaultaneoua  integral  equations  trust  be 
used.  The  purpose  of  this  Appendix  is  to  present  a  technique  for  dealing  with  a 
"squirter"  for  which  the  ratio  2h/a  is  not  small  enough  to  permit  setting  the  co¬ 
efficients  B  in  Eq.  EC. 10  equal  to  unity,  even  though  the  corresponding  ratio 
2h/X  is  sufficiently  small  to  allow  us  to  ignore  the  compressibility  of  the  fluid 
in  the  annular  region.  The  most  straightforward  approach  is  to  evaluate  the  un¬ 
perturbed  potentials,  which  are  obtained  by  setting  oc(z)  =»  0,  and  to  use  these 
potentials  to  compute  a  perturbation  solution  of  a(z)  from  Eq.  EC. 2: 

L 

a(°)(z)=  -  SSj  18o(1+  f)g0(atfao,2-z')  +  8^1-f  k^.o^z-z')]  dz'  (B.l) 

The  perturbation  solution  of  the  impedance  correction  factor  is  obtained 
by  substituting  this  expression  for  Ct(z)  in  Eq.  VHI.8,  with  T^z-z')  as  defined 
in  Eq-  EC. 14.  Tbs  far  field  potentials  can  of  course  also  be  obtained  in  »  straight¬ 
forward  fashion  by  substituting  ot^(z)  in  Eq.  VU.8  and  9>  for  4q,  and  Bq.  VECI.3 
tor 

'Die  perturbation  solution  can  be  improved  by  iteration  03  follows:  One  sub¬ 
stitutes  cr^°'(zO  for  a(z')  in  the  z'- integral  in  Eq.  EC. 11  and  solves  for  a(z). 

This  8 counts  tc  solving  Eq.  EC.l  for  ar(z)  using  the  perturbation  solutions  of  iQ 
and  fj..  If  this  Iteration  process  is  repeated  p  times,  one  finds  that  the  pth 


iterate  of  a(z)  is  related  to  the  (p-l)th  iterate  aa  follows: 


u'p;(z)  =  -  {J  [80(l+  +  e1(1-  ,*-»')]  dz' 

0 

CD 

+  J  a(p*l)(z/)  r(z-z')  dz'}  (B.2) 

where  r(z-z')  is  defined  in  Eq.  JX.lte. 

We  will  now  present  a  finite-difference  procedure  for  solving  the  integral 
equations.  Instead  of  requiring  that  the  integral  equation,  Eq.  IX. 11,  he  satis¬ 
fied  for  all  valusa  of  z  in  the  region  |z|  >  1,  ve  satisfy  it  at  a  finite  number 

of  points  zn  »  z^jZgjZ^ . Zj|  separated  by  intervals  2dQ.  These  intervals  should 

be  selected  saaller  in  regions  close  to  the  transducer  extremities ,  which  make  a 
more  important  contribution  to  the  potentials  than  more  distant  regions.  Further¬ 
more,  we  assume  that  the  unknown  function  a(z)  has  a  constant  value  oq  in  each 
interval  (z^-d^)  <  z  <  (z^+d)  and  varies  dis continuously  from  one  interval  to  the 
next.  We  thus  arrive  at  a  set  of  K  simultaneous  equations  in  H  unknown 
quantities  c  : 


2dir(o)  ^2r(v22>-  •  •  ®V<vV  M 


ai^Z2-Zl)  +  2d2r  <0) 


2dir(vsi) 


A  +  2V<°) 

a 


L 

F(sn}  =  J [0O(1+  J)60(Vb0»V8,)  +  9i(l*  427  (B*3) 

The  two  Green's  functions  which  enter  into  the  linear  const  lent  ion  r(z-z'),  Eq.  IX.l*ia. 
have  a  pole  at  z»z'.  The  diagonal  terms  in  the  above  matrix  do  not,  however,  display 
a  singularity,  3lnce  they  are  equivalent  to  an  Integral  of  r(z-z')  over  z',  vb5ch, 
like  the  potential,  is  well  behaved. 


Solving  this  set  of 
impedance  component 

Za  =  -2(lifca)2 


equations  for  the  values 
Za  from  Eq.  VHI.8 


of  aQ  we  can  compute  the  radiation 

(B.4) 


Robey's  impedance  cosponent,  Z^,  is  given  in  Sq.  jDC.19. 

When  applying  the  finite-dif  fereiso  aethod  to  the  non-axloytroetric  velocity 
distributions  discussed  in  Section  X,  it  is  not  necessary  to  construct  a  two- 
dimensional  grid  of  points  (zn,<pp)  over  the  tvo  cylindrical  surfaces  r=a, jz  |  >  L. 
Rather,  a  one-dleenslonal  set  of  finite-difference  equations  in  a,  of  the  fora  of 
Eq.  B.3  applies  to  each  iwdal  velocity  distribution  &Q(z,'  associated  with  the 
non-axisyiaasi  ic  Green's  functions,  Eq.  X.3  and  5. 

In  •«  v;,  it  is  recalled  (end  of  Section  IX)  that  a  variational  solution 
is  applicable,  oven  when  2h/a  is  not  but  that  Za  .-annot  be  obtained  direct¬ 

ly  froo  the  functional  J(a],  without  also  solving  for  the  unknown  coefficients  in 
the  trial  function  a(z). 


4l 


(Baron,  Matthews  and  Bleich,8 
Chen  and  Schwelhert3) 

Grid  of  point-sources  approximates 
radiating  surface.  Source  strength 
deter&ined  free  finite-difference 
solution  of  Helcholt2  integral  equation 


Robey's  Matheaatical  Model  of 


the 


(l*iird  and  Cohen,6  Robey,6 
Greenspen  ,7,s  and  ShernBn* ) 

Integral  equation  circuaveuted 
by  assuming  rigid  cylindrical 
baffles  r=a,  |  z  |  >  L,  and  by  con¬ 
structing  Green* 8  function  for 
which  (dG/dr')=0  for  r,=Ji 

l)  Known  velocity  distribution 
over  radiating  surface 

' ?)  Seed- infinite  rigid  baffles 


(c)  Robey's  Hatheaatical  Model  of  "Squirter" 


(Robey”  ) 

Seek  as  Pig.  l(t)  but  p^jane  baffles  in 
regions  r  >  a,  !  »  t  l 

iT)  Infinite  rigid  baffles 


Pig.  1.  REVIEW  OF  PUBLISHED  AEALTSES  OP  CILEtDRICAL  RADIATORS 

(See  Table  2  on  p.  3) 


mum 


Fig.  2.  SWttHX  OP  PK3S33T  APPROACH  (See  Table  2  on  p.  3) 


1)  Knovn  velocity  distribution  of  the  radiating  surface  (-L  <  z  <  L,  r=a) 

2)  Unknown  velocity  distribution  u(z)  satisfies  Integral  equation  on  sur¬ 
faces  r*a,  | s |  >  h  (z-dependent  pl»8e  shift  of  a(s)  is  not  indicated] 

Radiation  iapedaacc  o  +  ZQ 

where  -1  inpedance  computed  iron  Robey's  asthesatlcal  model,  Fig.  lb 
ZQ  ■  correction  associated  with  unknown  velocity  distribution  a(z) 
Variational  principle:  for  correct  a(z),  ZQ  <*  J(a] 

Ml2lt0 

fa 


CONTOUR  IIIT3G  RATION  OP  THE  INVERSE  TRANSFORM  OF  THE  GREEN'S  FUNCTION  G. 
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mvE  Underwater  Explosion  Res.  Div. 

Bor  folk  Kaval  Shipyard 
Portsmouth,  Virginia 
Attn:  Mr.  D.  S.  Cohen 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  335,  Tech.  Info.  Div. 

Code  31'5 ,  Mr.  F.  Vane 

Code  teO 

Code  tel 

Code  te3 

Code  te5 

Code  Vtf) 

Code  442 
Code  443 

Code  689,  Mr.  I.  Cook 
Code  1500 

Chief,  Bureau  of  Kaval  Weapons 
Department  of  the  Kavy 
Washington  25,  D,  C. 

Attn:  RAAD,  Airframe  Design 
DLI-3,  Tech.  Library 
R-12,  Chief  Scientist 
I*CA,  M3  and  Airframe  Br. 

HI,  ASW  Division 
RKHE,  Research  Br. 

Special  Projects  Office 
Bureau  of  Kaval  Weapons 
Department  of  the  Havy 
Washington  25,  D.  C. 

Attn:  Code  SP-001,  Chief  Scientist 
Code  SP-20,  Tech.  Director 

Chief,  Bureau  of  Yards  and  Docks 
Department  of  the  Havy 
Washington  25,  D.  C. 

Attn:  Code  70,  Research 

Code  E228  Tech.  Library 

Cosaanding  Officer  and  Director 
David  Taylor  Model  Basin 
Washington,  D.  C.  20007 
Attn:  Code  108,  Mr.  R.  T.  McGoldrick 
Code  10&,  Dr.  M.  Strasberg 
Code  140,  Tech.  Info.  Div. 

Code  538,  Dr.  F.  Theilhclmer 
Code  563,  Mr.  A.  0.  Sykes 
Code  700.  Dr-  A.  ”.  Kell 
Code  720,  Mr.  E.  E.  Johnson 
Code  731,  Kr.  J.  G.  Pulos 
Code  740,  Dr.  W.  J.  Sette 
Cede  760,  Mr.  E.  Hoonan 
Code  761,  Dr.  E.  Buchsan 
Code  771,  Dr.  R.  Liebovltz 

Commanding  Officer  and  Director 
U.S.  Kavy  Underwater  Sound  Laboratory 
Fort  Trumbull 
Hew  London,  Connecticut 


Cocnender 

U.8.  Haval  Ordnance  Laboratory 
White  Oak,  Maryland 
Attn:  HL,  Technical  library 
D,  Technical  Director 
<XI,  Uudervater  Weapons 
RS,  Acoustics  Division 

Cosaanding  Officer 

U.S.  Kaval  Kins  Defense  Laboratory 

Panama  City,  Florida 

Coacander 

U.S.  Saval  Air  Development  Center 
Johnsvllle,  Pennsylvania 

Director 

U.S.  Ravy  Underwater  Sound  Ref.  lab. 
Office  of  Haval  Research 
P.  0.  Box  8337 
Orlando,  Florida 

CoEnanding  officer  and  Director 

U.S.  Havy  Electronics  laboratory 

San  Diego  52,  California 

Attn:  Mr.  George  S.  Colcaon,  Code  2323 

Cocrander 

Portsmouth  Saval  Shipyard 
Port'.nouth ,  Hew  Hampshire 

Cosed nder 

te-oe  Island  Saval  Shipyard 
Vallejo,  California 

Director,  Materia  la  Laboratory 
Hev  York  Saval  Shipyard 
Brooklyn.  Hev  York  11251 

Officer- in-Charge 
Haval  Civil  Engineering  Research 
and  Evaluation  Laboratory 
U.S.  Haval  Construction  Battalion  Center 
Port  Hueneae,  California 

Director 

Haval  Air  Esperiaent  Static« 

Haval  Air  Material  Center 

Haval  Base,  Philadelphia  12,  Penna. 

Attn:  Structures  Laboratory 

Officer-in-Cterge 

David  Taylor  Model  Basin 

Underwater  Explosion  Research  Division 

Bor  folk  Haval  Shipyard 

Portsnouth,  Virginia 

Attn:  Dr.  H.  M.  Schauer 


Coracander 

U.  S.  Saval  Proving  Ground 
Dahlgren,  Virginia 

Supervisor  of  Shipbuilding,  U3H  and 
Haval  Inspector  Ordrance 
General  Dynamics  Corporation 
Electric  Eoot  Division 
Groton,  Connuctittit  063W) 

Coanander 

Haval  Ordnance  Test  Station ,  Cnine  Lake , 
Attn:  Phyeics  Division  California 
Mechanics  Division 

Ccweading  Officer 
Saval  Ordnance  Te3t  Station 
Underwater  Ordnance  Dl  *~..on 
3202  B.  Foothill  Boulevard 
Pasadena  8,  California 

Commanding  Officer  and  Director 
U.S.  Haval  Engineering  Experiment  Sta, 
Annapolis,  Maryland 

Superintendent 

U.S.  Saval  Postgraduate  School 
Monterey,  California 

Coasaander 

Air  Material  Cosaaand 
Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 
Attn:  HCREX-B 

Structures  Division 

Coanander ,  VJADD 

Wright-Patterson  Air  Force  Base 
Ohio 

Attn:  WWRC 
WWJMDB 
VWHKDD 

Director  of  Intelligence 
Headquarters,  U.S.  Air  Force 
Washington  25,  D.  C. 

Attn:  P.  V,  Branch  (Air  Targets  Dlv.) 

Ccoaander 

Air  Force  Office  of  Scientific  Res. 
Washington  25,  D.  C. 

Attn:  Mechanics  Division 

U.S.  Atomic  Energy  Cesaaissioa 
Washington  25,  D.  C. 

Attn:  Director  of  Research 
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Director 

national  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Division  of  Mechanics 

national  Aeronautics  and  Space  Adm. 
1512  H.  Street,  N.  V. 

Washington  25,  D.  C. 

Attn:  Chief,  Div.  of  Research 
Information 

Director 

Rational  Aeronautics  and  Space  Adm 
Langley  Research  Center 
Tnngley  Field,  Virginia 
Attn:  Structures  Division 

Dr.  M.  L.  Baron 

Paul  Weidlinger,  Consulting  Engineer 
770  Lexington  Avenue 
Nev  York  21,  Hev  York 

Professor  H.  H.  Bleich 
Deportment  of  Civil  Engineering 
Columbia  University 
618  Kudd  Building 
Hev  York  27,  *ev  York 

Professor  B.  A.  Boley 
Department  of  Civil  Engineering 
Columbia  University 
6l8  Kudd  Building 
Nev  York  27,  Mev  York 

Professor  Nicholas  J .  Hoff 
Dept,  of  Aeronautical  Engineering 
Stanford  University 
Stanford,  California 

Dr.  C.  W.  Horton 
Defense  Research  laboratory 
University  of  Texas 
Austin  12,  Texas 

Dr.  M.  A.  Brull 
Engineering  Mechanics  Division 
University  of  Pennsylvania 
Philadelphia  U,  Pennsylvania 

Dr.  ?.  DIKaggio 

Dept,  of  Civil  Engineering 

Columbia  University 

6l8  Kudd  Building 

York  27,  Hev  York 

Professor  D.  C.  Druoker 
Division  of  Engineering 
Brown  University 
Providence  12,  Rhode  Island 


Dr.  Ira  Dyer 

Bolt  Beranek  end  Newman,  Inc. 

50  Moulton  Street 
Casbridge,  Massachusetts  02138 

Professor  A.  C.  Eringen 
Dept,  of  Aeronautical  Engineering 
Purdue  University 
Lafayette,  Indiana 

Dr.  Martin  Goldberg 
Special  Projects,  Applied  Mechanics 
Grusnsn  Aircraft  Engineering  Ccrp. 
Bethpage,  Long  Island,  Hev  York 

Professor  J.  H.  Goodier 
Department  .  f  Mechanical  Engineering 
Stanford  University 
Stanford,  California 

Dr.  Josh  E.  Greenspon 
J  G  Engineering  Research  Assoc. 

3 — ,9  Callaway  Avenue 
Bnlti&ire  15,  Maryland 

Professor  Philip  G.  Hodge 
Armour  Research  Foundation 
10  West  35th  Street 
Chicago,  Illinois 

Professor  W.  Prager,  Chairman 
Physical  Sciences  Council 
Brown  University 
Providence  12,  Rhode  Island 

Professor  Joseph  Kempner 
Dept,  of  Aeronautical  Engineering 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
333  Jay  Street 
Brooklyn  1,  Hev  fork 

Professor  J .  M.  Klosner 
Dept,  of  Aeronautical  Engineering 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
333  Joy  Street 
Brooklyn  l,  Hev  York 

Professor  E.  H.  Lei 

Brovn  University 

Division  of  Applied  Mathematics 

Providence  12,  Rhode  Island 

Professor  R.  D.  Kindi in 
Dept,  of  Civil  Engineering 
Columbia  University 

6l3  Kudd  Building 
Rev  York  27,  Hev  York 
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Professor  P.  M.  Kaghdi 
University  of  California 
College  of  Engineering 
Berkeley,  California  04-0371 

Professor  H.  M.  Hevmark,  Head 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois 

Professor  P.  Pohle 
Department  of  Mathematics 
Adelphi  College 
Garden  City,  Rev  York 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  !&ssachusetts 

Professor  F.  V.  Ruaano 
Dept,  of  Aeronautical  Engineering 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
333  Jay  Street 
Brooklyn  1,  Rev  York 

Professor  E.  Reiss 

Institute  of  Mathsaatical  Sciences 

Hew  York  University 

25  Waverly  Piece 

Hev  York  3,  Rev  York 

Dean  V.  L.  Salerno 
College  of  Science  and  Engineering 
F*«irleigh  Dickinson  University 
Tea neck,  Hev  Jersey 

Professor  A.  S.  Velatso3 
Department  of  Civil  Engineering 
University  of  Illinois 
Urfcana,  Illinois 

Brovu  University 
Research  Analysis  Group 
Providence,  Rhode  Is lend 

Hydrospace  Research  Corporation 
17*9  Rockville  Pike 
Rockville,  Maryland 

Hudson  Laboratories 
Columbia  University 
1*5  Palisades  Street 
Dobbs  Ferry,  Hev  York 

Lament  Geological  Observatory 
Columbia  University 
Torre  Cliffs 
Palisades,  Hev  York 


Ordnance  Research  Laboratory 

University  of  Pennsylvania 

P.  0.  Box  30,  State  College,  Pa.  16601 

Dr.  Richard  Waterhouse 
Aoericnn  University 
Physics  Department 
Washington  16,  D.  C. 

Mr.  Stan  Lemon 

Chesapeake  Instrument  Corporation 
Shadyside,  Maryland 

Dr.  M.  M.  Backus 
Texas  Instruments  Incorporated 
100  Exchange  Park  Kortb 
Dallas,  Texas 

Dr.  H.  Basin 
Hughes  Aircraft  Company 
P.  0.  Box  2097 
Fullerton,  California 

Mr.  John  Mahoney 

OHS  Resident  Representative 

Harvard  University 

*73  Broadway 

Cambridge,  Mass.  02138 

Kaval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland 
nttn:  Dr.  D.  ?.  Bleil 
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